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In the light of experience gained over the course of the past year it is now possible 
to make a clearer statement of the basis cn which division into the two Sections. 
of the Proceedings is being made. 

In general, it has been found that the fektar which best meets the needs of 
Fellows is one which brings papers on microphysics and the physics of elementary 
particles into the one Section, A, and papers on macroscopic physics into the 
other Section, B. While it is essential to maintain a certain amount of flexibility 
in the allocation of papers, a paper will, as far as possible, be put in the Section 
where the main interest lies, e.g. a paper on counter-technique would appear in 
B if the circuitry is the main interest, but would be in A if the applications to actual 


counting are of major importance. 


A list indicating this subject division is given below. 


A.C. STICKLAND, 
Secretary—Editor. 
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Colour physics 

Astrophysics and solar physics 

Applications of spectroscopy (e.g. 
analysis, structure analysis) 

Analysis of crystal structure 
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The Elastic Constants of a Solid containing Spherical Holes 


By J. K. MACKENZIE * 
H. H. Wills Physical Laboratory, University of Bristol 


MS. received 4th August 1949 


ABSTRACT. The effective bulk and shear moduli are calculated by a self-consistent 
method due to Fréhlich and Sack. The bulk modulus k is determined by applying a 
hydrostatic pressure, and the shear modulus p by applying a simple homogeneous shear | 
stress, to a large sphere. Each hole is surrounded by a spherical shell of real material, and 
the reaction of the rest of the material is estimated by replacing it by equivalent homogeneous _ 
material. For consistency, both the density and the displacement of the outer spherical - 
boundary must be the same whether the hole and its surrounding shell are replaced by 
equivalent material or not. The effective elastic constants calculated from these conditions 
are 
1/k=1/Rop + 3(1 —p)/4pop + O[(1 —p)*], 

(Ho — #)/Mo= 5(1 — p)(3Ro +4 p40)/(9o + 840) + OL — p)*I, 
where ky and py refer to the real material and p is the density of the actual material relative 
to that of the real material ; in the next approximation k depends on the standard deviation 
of the volumes of the holes. 

The dilatation due to a distribution of pressures in the holes is (1/k—1/ky), where p 
is the mean obtained when the pressure in each hole has a weight proportional to the volume 
of the hole. By using the hydrodynamic analogue of the elastic problem, the theory is 
briefly applied to the theory of sintering, and used to discuss the effective viscosity of a 
liquid containing small air bubbles. 


§1.. INTRODUCTION 

‘HILE developing a theory of sintering (Mackenzie and Shuttleworth 

1949) it became necessary to calculate the compressibility of a solid 

containing spherical holes. In this paper all the holes were assumed 
to be of the same size and embedded in an incompressible medium. In this 
simple case it suffices to equate the stored elastic energy to the work done by the 
externally applied pressure in order to derive the effective compressibility. The 
original calculation, however, followed the lines of ordinary elasticity theory, 
and both the effective compressibility and the effective shear modulus were 
derived without assuming either holes of equal size or incompressibility of the real 
material. It is these calculations which are reported in the present paper. On 
the other hand, the present method can only be used when the relation between 
stress and stain is linear, while the energy method used in the paper on sintering 
is applicable when this relation is non-linear. 

In §2 the problem is analysed in some detail and the self-consistent method 
of calculation is described. In §3 two problems are considered: (a) the calcula- 
tion of the compressibility of a material full of small spherical holes of different 
sizes, and (bd) the calculation of the dilatation due to different pressures acting in 
the holes; the shear modulus of a material full of small spherical holes is calculated 
in§4. Finally, in §5 the close similarity between the equations of hydrodynamics 
and elasticity is used to discuss briefly an application to the theory of sintering and to 
discuss the effective viscosity of a liquid containing small air bubbles. 

The results of these calculations show that, to a good approximation, both the 
effective elastic moduli depend only on the relative density p and the elastic 


* Now at C.S.I.R.O., Division of Tribophysics, Melbourne, Australia. 
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constants of the real material; holes have a comparatively large effect on com- 
pressibility and only a small effect on the shear modulus. It may be conjectured 
that this is true in general no matter what the shape or size of the individual holes 
may be, provided of course that extremes are excluded. It is not clear how the 
numerical constants appearing in the results will depend on the particular shape 
of hole chosen for calculation, but it seems likely that the formulae obtained here 
will give the correct order of magnitude in most cases. 


Ses) DH PME THOD TOR SCAT CU MAT ON 

It will be assumed that the real material contains isolated spherical holes 
distributed at random throughout the volume of the material, and that the real 
material has homogeneous and isotropic elastic properties. It will also be assumed 
that the volume of all the holes is small compared with the total volume, and that 
the total volume contains a large number of holes. The calculation uses the notion 
of an equivalent homogeneous continuum introduced in the next paragraph. 

If attention is directed not to the detailed structure of the actual solid but to 
avolume element containing a large number of holes, then, because the distribution 
of the size of holes and of their position in space is statistically uniform, such a 
volume element can be regarded as forming part of a homogeneous isotropic 
continuum. The elastic constants of this equivalent homogeneous continuum 
will be calculated in terms of the elastic properties of the real material and the 
number and size distribution of the holes. When large volumes are concerned, 
the elastic properties of the actual solid will be independent of its external shape, 
and for simplicity it is convenient to choose this shape to be spherical. The 
effective elastic constants will be calculated by applying a stress to the outer 
spherical boundary and comparing the strains produced in the actual solid with 
the strains produced by the same stress in a homogeneous isotropic solid. A 
uniform hydrostatic pressure will be applied to determine the effective bulk 
modulus and a simple homogeneous shear stress (with no hydrostatic component) 
to calculate the effective shear modulus. The results obtained depend to some 
extent on the particular stress whichis applied and, in general, the more complicated 
the stress the more the detailed structure of the actual solid will become apparent. 

A convenient self-consistent method of calculation has been indicated by 
Frohlich and Sack (1946). Their method is equivalent to a perturbation calcula- 
tion, and is certainly valid when the volume of all the holes is sufficiently small 
compared with the total volume. ‘The method depends on the fact that any 
macroscopic volume (i.e. one containing a large number of holes) can be replaced 
by equivalent homogeneous material without essentially changing the elastic 
behaviour of the solid. This is so because the mean stresses and displacements 
at the boundary of the volume are equal to those at the boundary of the same 
volume in the equivalent elastic continuum. The conditions for consistency, 
namely that the density and the displacement of the outer boundary shall be 
unchanged by such replacement, enable the effective elastic constants to be 
calculated. The particular method of calculation adopted is described in detail 
an the next paragraph. 

Since the holes are isolated, every hole of radius a will be surrounded by a 
spherical shell of real material out to at least some radius 7) (see Figure), and this 
may be regarded as a typical volume element of the actual solid. ‘The model used 
for the calculation of the effective elastic constants is then as follows. A hole of 
radius a is surrounded by aspherical shell of real material out to radius 7y, and this 

A-2 
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in turn is surrounded by a spherical shell of equivalent homogeneous material» 
(whose elastic constants are as yet unknown) out to somelarge radius R. A particular 
stress is applied to this outer spherical boundary and the displacements of the outer 
surface calculated ; this displacement is then weighted according to the frequency — 
of occurrence of holes of radius a and averaged over the size distribution of the 
holes. Finally, the resulting mean displacement is equated to the displacement | 
which would have resulted had the sphere of radius R been completely filled with | 
equivalent homogeneous material. Thus, one of the consistency conditions is | 
satisfied and an equation for the effective elastic constants derived. 


ae 


Ma 


eal- 
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Equivalent Homogeneous 
Continuum 


The model used for calculation. 


The other condition determines 79, since the density of a hole and its surround- 
ing shell must be the same as that of the actual material. If f(a)daisthe proportion 
of holes with radii lying in the range (a, a+da) then the volume of holes per unit 
volume of the real material is 


— = an | faa’ da= | faa’ da; =~ “eee (1) 


where the integration is over the whole of the relevant range of a, p is the relative 
density, and 7 is the number of holes per unit volume of the real material. Note 


that | fla) da=1. 


In general, when part of a solid is replaced by material of different elastic 
constants there is an additional term in the displacement of the outer boundary 
of the order of the ratio of the volume replaced to the total volume. If, however, 
the elastic constants of the new material are the same as the effective constants of 
the original material this term will disappear. In fact, the consistency condition 
will be expressed by equating this term of order 793/R® to zero. In addition, there 
will usually be terms of higher order which arise because the boundary conditions 
are never met exactly when part of the actual solid is replaced by equivalent 
continuum, but only on the average. ‘These terms will be neglected. 

In the following sections the method outlined above will be applied to the 
calculation of the bulk and shear moduli of a material containing small spherical 
holes. 
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$3. THE COMPRESSIBILITY 

Two problems will be considered in this section: (a) the calculation of the 
compressibility of a material containing small spherical holes of different sizes, 
and (6) the calculation of the dilatation due to a distribution of pressures inside the 
holes; these pressures may be due to included gas or to surface tension in the 
boundary of the hole. The compressibility will be calculated by applying a 
hydrostatic pressure P to the outer spherical boundary and using the method 
described in the last section. ‘The same model is also applicable when a pressure 
p is applied inside the central hole with zero pressure in all the other holes ; in this 
case, the approximation consists in estimating the reaction of the part of the actual 
solid which lies outside a sphere of radius ry by replacing it by equivalent homo- 
geneous material. When there are pressures in all the holes the total displacement 
is found by adding the displacements due to each hole individually. 

The detailed calculations fall naturally into three parts. In §3 (1) the general 
equations for the elastic displacements will be written down, while in the two 


_ following sub-sections the effect of an externally and an internally applied pressure 


will be considered in turn. However, before proceeding it is convenient to show 
how the second problem may be reduced to the first. 

Consider first the case where the pressure p is the same in every hole and there 
is no externally applied pressure. If an additional hydrostatic pressure —p is 
applied to the surfaces of all the holes and to the external surface, then the resultant 
pressure is zero in all the holes and —p on the external surface. Now when a 
hydrostatic pressure is applied to all the boundaries of a homogeneous elastic 
material the stress at all points inside the material is a hydrostatic pressure equal 
to that applied (Love 1944, p.86). Hence, the displacements due to the extra 
hydrostatic pressure —p can readily be estimated by filling up all the holes with 
real material and applying the pressure —ptothe outer boundary. (A comparison 
of the two situations shows that all the boundary conditions are met.) ‘Thus the 
total dilatation is p(1/k—1/k,), and the additional term is zero if the material is 
incompressible. 

On the other hand, if the pressure is not the same in all the holes, then a fixed 
part of this pressure can be transferred to the outer boundary as explained in the 
last paragraph, and the effect of the pressures remaining in the holes estimated 
separately. If these remaining pressures produce no total dilatation they may be 
neglected. Now the displacement of the outer boundary when only one particular 
hole has a pressure in it is proportional to the pressure in the hole and to its volume. 


_ Therefore, if the pressure in each hole is yiven a weight proportional to the volume 


of the hole, then when the mean pressure so calculated is transferred to the outside 
the remaining distribution of pressures will produce no displacement of the 
external boundary. ‘This statement will be verified in § 3 (iii). 


(1) General Elastic Equations 


First a word about notation. Symbols with subscript zero refer to the real 
material or to the intermediate boundary, while the same symbol without a sub- 
script usually denotes the corresponding magnitude for the equivalent homo- 
geneous continuum. Thus A and vw are the effective elastic constants which are 
to be calculated, while X, and yz, are the same constants for the real material; 
the constant yu is the rigidity and the bulk modulus k=A + 3p. 
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It is required to find the displacements in the composite body illustrated in 
the Figure when hydrostatic pressures are applied to both boundaries. On 
account of the spherical symmetry, the displacements will be purely radial and | 
given by (Love 1944, p. 142) | 

u=Ar+B/r’, | 
while the normal pressure acting across the spherical boundary at the same radius Is. 
Ou u | 

(A+ 2p) = +2A-. 
Then, using Love’s results for a thick spherical shell, the radial displacement in 
the outer shell of equivalent homogeneous material is =| 


Opie ee Ry? py — P 
“3k SR 78 4ur* R—7,3’ 


u 


and the radial displacement in the inner shell of real material is 


eR S pote Glee P ae 
= ape TS , 
Ske fo oe Auer? 1° — a? 


Uu 


where p is the pressure inside the hole (all other holes at zero pressure) and P is 
the externally applied pressure. 

The pressure p, at the intermediate boundary is determined by the boundary 
condition that requires the radial displacement to be continuous across this 
boundary. This condition gives, neglecting 793/R? and higher powers, 


Ly Side Whew. ae Ee eee eee tte latalece 
Pol 4u.? 3h, 78a? day Fo Sat) | OB aa le eel 


This equation for py, together with the expression for the displacement of the 
external boundary, namely 


u 


P oth ate 
te EGP) ee 3 
R= 3p {Po ae ae ©) 


completely determines the external change of volume. The detailed solution of 
the two problems can now be obtained quite simply. 


(ii) Pressure applied Externally 


In this sub-section the effective compressibility k is determined when there 
are no pressures inthe holes. Whenapressure P is applied to the external surface, 
the resulting displacement of the external boundary due to holes of some particular 
radius a is given by equation (3). Averaging this displacement over the size 
distribution of the holes and equating the resulting mean displacement to the 
displacement of the external boundary when a pressure P is applied to a solid 
sphere of the equivalent homogeneous material, viz. — PR/3k, as required by the 
consistency condition, it follows that 


P= | p)f(a)da, = eens (4) 
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where fp is given by equation (2) and | Aa) da=1. If the holes are all the same 
size, py =P, and (7) (with o=0) follows from (2) immediately. 
On putting p=0 into equation (2) and expanding in powers of a3/7,°, 


Po __ 1+4u/3k @& feels S =} + | 


P= (4u3kl ore 1 4p 3k, 
Now it will be shown in §4 that 


Gi Pea p tb ange Tateiers (5) 
where « depends only on ky and py and is given by equation (20). Put 


{Sha ate eee eee eth Se ot eee (1) 
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and 
[Fahlaydamot+(l—p see (6) 


where 4707,°/3 is the standard deviation of the volumes of the holes. Then it 
follows trom (4), after some reduction, that 


This equation gives the effective compressibility up to and including terms of 
order (1—p). The value of x is always about 2. It should be noted that a few 
holes in an otherwise incompressible medium can give rise to an appreciable 
finite compressibility. 

It is not altogether certain that the model is adequate for the calculation of a 
third approximation to k. It is clear, however, that the effect of a distribution 
in the sizes of the holes is small, and that to a good approximation the effective 
compressibility depends only on p and the elastic constants of the real material. 
In fact, (7) gives the correct behaviour for the effective compressibility in the 
extreme cases as the density tends to unity or to zero. Thus it is plausible that 
this equation should be roughly correct for all relative densities. 


(iii) Pressure applied Internally 

In this sub-section the effect of a distribution of pressures inside the holes. 
is considered. If there are different pressures in holes of the same size it is clear 
that their effect is the same as the mean pressure p acting in all of them. | When 
a pressure p is applied in only one hole, equation (3), with P =0, gives the resulting 
displacement of the external boundary. Then, on multiplying by the number 
of holes, R?/r,?, and averaging over the size distribution, the mean displacement 
of the external surface is (neglecting r,?/R® compared with unity) 


u=R E ze ral | Pof(a)da=pR E = | a eee (8) 


where p is the mean pressure which is transferred to the external surface, and the 
last part of the equation follows from the discussion at the end of the introductory 
part of §3. ‘The last equation in (8) determines the mean effective pressure 
p, which will now be calculated. 
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On putting P=0 into equation (2), substituting the value of py into (8), and 
expanding in powers of a3/7,%, it follows that 


~_ 1/3ky + 1/414 ite | =| 14 (Ho — #)/Ho ©; |fayaa. 


~ -1/3k—1/3k, 14+4/3k, 1+41/3Ry 76° 


‘Then, substituting for k from (7) and replacing p by p when evaluating the 
correction term in the integral, it follows, after some reduction, that 


p= a | poi ada 4 ie a ene (9) 


where a,° =(1 —p)ro?, or 477a,°/3 is the mean volume of the holes. Although this 
expression is correct up to and including terms of order (1—p)%, it is not clear — 
what error has been committed in replacing p by p in the small correction term. 
Equation (9) justifies the statement made at the very end of the introductory part 
of $3. 

Finally, returning to equation (8), it is clear that the dilatation caused by a 
distribution of pressures in the holes is p/k’, where p is given by (9) and 


P7a-R7 oe le tet oe  eiteke anette ( ) 


§4. THE SHEAR MODULUS 

The effective shear modulus will now be calculated by applying a simple type 
of homogeneous shear stress and calculating the resulting displacements of the 
external boundary both for the composite solid shown in the Figure and for a 
sphere of equivalent homogeneous material. As before, comparison of the two 
solutions gives an equation from which the effective shear modulus can be calculated. 

Now the displacements which result from all the simple types of homogeneous 
shearing stress can be derived from solid spherical harmonics of degree two. 
Such a solid harmonic is of the form r2.8,, where S, is a surface harmonic of degree 
two and determines the symmetry of the displacements; the solid harmonic 
r*S. also has the same symmetry, and will therefore be expected to occur in the 
general solution. 

The solution of the elastic equations, which is derived from the solid harmonic 
7°S., and for which the dilatation and the rotation are everywhere zero, is (Love 
1944) p25) 

Uperad(7*S,)=27S,6rS5-)) ee eee (11) 


and the corresponding stress across a spherical boundary at the same radius is 
Be Jur grad @2S85)=4uS—r + 2a5, vy emen e (12) 


where r and + are unit vectors normal and tangential to the spherical boundary 
and S,’ is independent of 7. This solution gives the displacements when a 
shearing stress F, which is the same at all points, is applied to a solid homogeneous 
sphere. This stress will now be applied to the composite solid and the resulting 
displacements calculated. 

There are four independent solutions of the elastic equations which can be 
derived from solid harmonics with the same symmetry as S,, and for which the 
rotation is everywhere zero. Writing 


pS.ta- GS n=) Gh = ee ee (13) 
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cthese solutions are (Love 1944, p. 250): 


u, = grad (r?S,) ={2, 1}r, 
u, =grad(r-3S,) ={—3, 1/74, 
6A 
u,=r'grad(r2S,)+o 73S. = oes i} a ee (14) 


WF ptad(?*S,) +047 oof = - aH ; i}/s ; 
fe 


the corresponding stresses acting across a spherical surface of radius r are, 
tespectively, 


F,/u={4, 2}, 
F,/u = (24, =8i/7°, 

SON ELE rll ry ae ora (15) 
Fule= {arpa SA+ 7p ee 


F,ju= i he, = J 


"Thus the general solution of the elastic equations in a thick spherical shell will 
be of the form 
G=]2077 bu, Cust DU. ae (16) 


‘This solution will apply as it stands in the outer shell containing equivalent 
homogeneous material; the corresponding solution for the inner shell containing 
the real material is derived simply by putting a subscript zero on all the constants. 
"The whole solution then involves eight arbitrary constants which are determined 
by eight equations arising from the boundary conditions; two equations have to 
be satisfied for each of four boundary conditions. 

The four boundary conditions are: (a) that the stress across the external boundary 
is u{4,2', (b) that the displacements are continuous across the intermediate 
boundary, (c) that the stress is continuous across the intermediate boundary, and 
(d) that the innermost boundary is stress-free. Using equations (14), (15) and 
(16), these boundary conditions lead to a set of equations whose formal solution 
can be written down immediately in a form involving eighth-order determinants. 
"The solution can then be expanded quite simply in powers of r/R, and the result 
-of doing this is 

(Pe 


ae eK: Re” Bik == 0; 
oes Reg a, © ears (17) 
Da I 8 10. = ma 0 3 
Ck ae R?? D A; R?? 


where powers of 7,/R higher than the third have been neglected, A, and A, are 
‘certain determinants of the sixth order, and K,, K, involve only A and pw. It is 
clear from this solution that the additional terms in the displacement due to the 
presence of the hole and its surrounding shell are of order r3/R*, as was stated at 


the end of §2. 
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Substituting the values of the constants into (16), integrating over the size 
distribution, and comparing the results with the displacement given by equation. 
(11), it follows that the condition for consistency is 


eae 
| xa) aa=0. Rr (18) 


Finally, expanding A, and A, in powers of a/r), and neglecting powers higher than. 
the third in the quotient A,/A,, it follows that 


io PE ay ey ee (19) 
Ko » Oh, + 8p 


This equation gives the fractional decrease in the modulus of rigidity due to the: 
presence of the holes. The value of the constant « introduced in equation (5) is. 


3h + 449 
rien ee re 


§5. APPLICATIONS 

The application of the results of the preceding calculations to determining the 
effective elastic constants of materials containing a large number of small holes 
is obvious and will not be discussed further. There is, however, a close similarity 
between the equations of elasticity and those of hydrodynamics (Goodier 1936), 
and it is the results which can be deduced from this similarity which will be 
briefly discussed in the following paragraphs. 

If, in the hydrodynamic equations for slow steady motion (so that D/Dt=0), 
the velocity v is replaced by the displacement u, the viscosity 7 by the rigidity p,. 
and —p by kdiv u, then the equations of elasticity are obtained, and vice versa. 
Thus, by means of the above scheme of replacement, the velocities in a slow and. 
steady viscous motion due to some system of stresses can be derived immediately 
from the elastic displacements resulting from the same stresses in the corres- 
ponding elastic problem. When making this transfer the bulk modulus will 
usually be taken as infinite because, although the shearing viscosity of a liquid 
is finite, the ‘volume viscosity’ (which is the analogue of k) is usually considered 
to be infinite, i.e. there is no energy dissipated during a pure compression or 
expansion. 

In the theory of sintering mentioned in the introduction, it was necessary to 
calculate the rate of increase of density of a compact consisting of spherical holes. 
in an incompressible medium when a negative pressure —2y/a due to surface: 
tension acts in al] the holes. In the case of glass the action of these pressures is. 
resisted by viscous forces, and when this is so the resulting flow may be derived 
from the results of §3(iii) by means of the above substitution. Thus the rate of 
increase of density is 


dp __pp 
Fre eee Mi oe (21) 
where p is given by equation (9) as 
INS Fie 
p= - 24 | ef(a)da | 
OE Bey th be? AD SG ye (22), 


2 mean surface area of holes 
3” mean volume of holes 


| 


a ae 
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If all the holes are of the same size, then using (10) with ky = 00 to determine k’,. 


dp 3 (42\* yn ‘ 
Ss = (5) cay =p ee > ee (23) 
a formula previously derived in another way; n is the number of holes per unit 
volume of the real material. 

Again, consider a liquid which contains a large number of small bubbles 
containing gas. These bubbles will be of such a size that the gas pressure is 
equal to the hydrostatic pressure in the liquid plus the pressure 2y/a due to surface: 
tension. Now provided the liquid is sheared slowly and the bubbles retain 
approximately their spherical form, the analogue of equation (19) gives the fractional. 
decrease in the viscosity due to the presence of the bubbles as 


i id big, OS Se aed (24) 
10 3 

where « is the fraction of the total volume which is occupied by the bubbles.. 
This expression assumes that the bubbles deform freely during the motion so that,. 
unless the motion is oscillatory, the bubbles will not remain spherical for long. 
In practice, both the gas pressure and the surface tension will tend to maintain 
the shape of the bubble approximately spherical and equation (24) will not apply. 
Taylor (1932) has shown that the presence of liquid spheres which remain 
spherical during the motion raises the effective viscosity rather than lowers it.. 
If the viscosity of the liquid in the spheres is y’, Taylor shows that 


To=79 _ 3 1 +2m/5 

No 2 9° +% 
When the spheres are rigid, 7’ = 00, and this reduces to Einstein’s result (1906,. 
1911), and when 7’ =0, 


a result also obtained by Eisenschitz (1933) for the case where the liquid slips. 
freely over the surface of a rigid sphere. It is clear therefore that boundary 
conditions assumed at the surface of the bubble play an important part in the: 
results which are obtained. 
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-ABSTRACT. Ferromagnetic particles of diameter from 10~* to 107° cm. when exposed 
to light move along the lines of force of an inhomogeneous magnetic field. It does not seem 
~possible to explain the movement when only dipole properties of the particles are taken into 
account. Field strengths of 10~? gauss are sufficient to cause distinct movement. ‘The 
~motions resemble those of electrically charged particles in electric fields. 


STS IN LRO DUCTION 

MALL particles of iron powder show, when exposed to light in homogeneous 
magnetic fields, polar motions which are reversible with reversal of the field 
(magnetophotophoresis (Ehrenhaft 1930)). Since then many kinds of 
various substances have been found to show this behaviour. In continua- 
tion of researches on the phenomenon the following observations were carried 
out in simply formed, typical inhomogeneous fields. [he method described for 

producing the particles is one of the most convenient and effective. 


§2. EXPERIMENT 
An electric arc between two iron electrodes in oxygen or air produces a sort 
-of smoke consisting mainly of small particles of ferric oxide. Such particles of 
10-* to 10-° cm. in diameter were collected in a small observation chamber which 
was surrounded by a small rectangular loop (3 mm. x 10 mm.) of thin copper wire 
(see Figure 1). ‘The chamber was specially constructed to allow of illumination 
-and observation near the wire. The particles were observed in a dark field by 
means of a microscope (magnification 30 to 100 times), the axis of which was 
parallel to the plane of the above-mentioned loop. The assembly consisting of 
-chamber, microscope and surroundings was iron-free and the geomagnetic field 
was compensated over the region of observation. 


§3. OBSERVATIONS 

Under the conditions of the experiment most of the particles appear as diffrac- 
tion discs, others have the form of short, often crooked, threads. They all fall uni- 
formly in the resisting medium (air) with a Brownian movement superimposed. 
In very few cases is common photophoresis (positive or negative) visible. 

As soon as a direct current flows through the wire of the loop (for instance 
0-2 amp.), forming an inhomogeneous magnetic field in the cell, 10 to 20°% of the 
particles visible as discs start moving along the magnetic lines of force, some in, 
some against, the direction of the magnetic field. They pass the plane of the 
loop unhindered. When the direction of the current is reversed, these particles 
instantaneously reverse their direction of motion. They stop immediately the 
current is cut off. ‘The other particles, so far as they appear as discs, do not show 
any reaction to the weak magnetic field applied. 

All the particles shaped like short threads behave like tiny compass needles: 
pointing in the direction of the acting magnetic field they turn through 180° 
when the field is reversed. Among these particles also, some show the above 

»described phenomenon, that is to say they move along the magnetic lines of force 


i 
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Figure 3. Time of exposure 10 sec. 
Magnification 40 times. 


Figure 4. Time of exposure 2 sec. 


Total magnification 120 times. 
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Figure 5. Time of exposure 0-4 sec. 
Magnification 80 times. Sunlight. 


Figure 6. Time of exposure 13 sec. 


Magnification 80 times. 


Movements of Particles in Magnetic Fields iis 


and pass the plane of the loop unhindered; reversal of the magnetic field causes 
a turn through 180° and reversal of the translatory movement. 

Taking into account only dipole properties, from a theoretical point of view 
the particles should move from the right and the left towards the plane of the loop: 
and the wire (provided that the gradient of the field is sufficient) without reversing 
their translatory movement. 

The apparatus as above described did not permit the use of higher values of | 
the inhomogeneous magnetic field. For this purpose the loop was replaced by a 
coil with iron core adjusted outside of the cell as shown in Figure 2, Small field 
strength (for example the remanent magnetism of the iron core) then caused 
movements analogous to those already described. 


GLOMAGN. FIELD 
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Figure 1. Figure 2. 


When now the magnetic field strength was gradually increased, the different” 
kinds of particles were affected in the following way : 


(i) All particles which showed no translatory movement at first began gradually 
(from about 20 gauss/cm.) to move towards the denser lines of force (i.e. towards 
the magnet). 


(ii) Those particles which were at first attracted by the magnet continued to 
move in this direction at an increasing velocity. 


(iii) Those particles which were at first repelled by the magnet gradually 
decreased their velocity, came to a standstill, and finally also moved towards the 
magnet with gradually increasing velocity. 


The same results can be achieved by steadily approaching a small permanent 
bar magnet to the cell instead of using a coil with iron core and increasing the- 
current as described above. 


§4. TRACKS OF PARTICLES 

Figure 3 is a general view showing tracks of some particles (the plane of the 
loop is vertical in the middle of the photograph). It is remarkable that the tracks 
often appear as dotted lines, which are in some cases enlarged to helix-like paths ; 
this is illustrated by Figure 4, which shows the track of a single particle, and by 
Figure 5. Examination by means of a stereo-microscope proved these motions 
to be essentially rotatory, and not a sort of fluttering due to translatory motion. 

Figure 6 shows for example the movements of a particle while the direction 
of the current in the loop was twice reversed. (Another particle showing only” 
common photophoresis crossed the tracks of the first-mentioned particle.) 
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§5. FURTHER OBSERVATIONS 


Ionization of the air in the cell by means of y-rays did not influence the effect 
in any visible manner. Common radiometer forces play apparently only a small 


part in these movements, because the effect of common photophoresis was found | 


to be only about 5 to 10% of the whole. 


With care the same particle could be observed over a long period (half an hour | 


and more), in course of which it did not lose the described property. Thus the 


explanation of the effect by assuming a jet-like emission of occluded gases caused _ 


by the impinging radiation, or as the result of a chemical action, cannot be accepted. 


{ 


Moreover, the velocity of the particles (in the direction of the magnetic field) is | 


not constant, but increases with increase of field strength (up to about 120 gauss). 


On an average the sensitivity of the particles in question is such that the— 
method can be used for determining direction and intensity of weak magnetic 


fields (e.g. the geomagnetic one) when an accuracy of only + 0-03 gauss is required. 


§6. CONCLUSION 


The described behaviour of the particles reminds one of equivalent experi- 
mental results, when electrically charged particles in inhomogeneous electric 
fields are observed, where the movements are described by the well-known 
equation 

Peek +flesk ye grad(Z*), 6 a ee (1) 
in which P is the force acting on the particle, e the electric charge, E the electric 
field strength, and fa function of the dielectric constant « and the shape K of the 
particle. (Conductivity is neglected.) 

Analogy suggests a trial of the equation 


P=GH +f(u, K). grad(H),. — © eee (2) 


where H is the magnetic field strength and « the magnetic permeability, but the 
‘coefficient of proportion q turns out to be a function of H, of intensity and compo- 
sition of light. In addition to this a remanent dipole momentum has to be taken 
into account in (2), when ferromagnetic substances are used. 

It is difficult to calculate values for g, because the shape of the particles (which 
can scarcely be assumed spherical as electron-optical photographs show) modifies 
Stokes’ Law ina rather uncalculable manner. The values of g are approximately 
10° (c.G.s.) and smaller, under the above-described circumstances. The effect 
is not restricted to irregularly shaped particles; it was found to be of the same 


order of magnitude with iron powder (precipitated from carbonyl iron), the 
particles of which consist of spheres. 
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ABSTRACT. Experiments on Geiger counters are described in which localization of the 
ion sheath was brought about by reducing the wire potential immediately after each count. 
"The effect on the dead time is discussed. The dead time of a counter 30 cm. long was 
found to be 20 usec. and independent of the counting rate up to 1-8 x 10° counts per minute. 
"The reduction of the dead time of short counters was also investigated. 


Si ENED INOD Ui Ciel ON 

T is an inherent property of the Geiger counter that each discharge is followed 

by an insensitive period or ‘dead time’, which may be of several hundred 

microseconds duration. This has been defined as “the time interval, after 
recording a count,....(during which) the counter is completely insensitive and 
does not detect other ionizing events occurring inside it”? (Korff 1946). Counts 
occurring during this period are missed and the registered counting rate must 
be corrected for their loss. The form of this correction is uncertain; Blackman 
and Michiels (1948) quote three different formulae for the correction factor and 
show theoretically that only one of them is valid under normal conditions. Further 
difficulties are introduced by the dependence of the dead time on the operating 
conditions, which may include the counting rate. 

One way of evading these difficulties is to work with such low counting rates 
that uncertainties in the correction to be applied may be neglected. Hence 
for 10° losses the counting rate must be limited to 10* counts per minute, if the 
paralysis time of the recording equipment is 600 psec. At higher counting rates 
the difficulty of applying an accurate correction increases rapidly. An empirical 
method of obtaining the true counting rate, when the dead time varies with it, 
has been described by Kohman (1945). This method is, however, rather laborious. 

Another expedient is the addition to the recording equipment of an electroni- 
cally determined paralysis time, much longer than the dead time of the counter. 
The purpose of the paralysis time it so reduce spurious counts and eliminate 
variations in the counter dead time. However, this does not solve the problem 
because the effective dead time is still not constant. This isso because asubsequent 
discharge of the counter, during the paralysis time, may prolong it by as much as 
the dead time. 

For many experiments the upper limit set to the counting rate by these 
considerations is too low. It would be desirable, therefore, to obtain an improve- 
ment, either by reducing the dead time, or by ensuring its constancy. bus, 
reducing the dead time would allow higher counting rates to be used for a given 
percentage of losses, while with a constant dead time larger corrections for losses 
could be safely applied. It would, of course, be most desirable to improve both 
factors. 

The dead time may be made constant by the use of a circuit which lowers the 
wire potential below the threshold of the Geiger region for a few hundred micro- 
seconds immediately after each count. The counter is then truly ‘ dead’ during 
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this constant period. There are several advantages of such a method of operation 


and Putman (1948) has shown that self-quenching counters operated in this way 
show improved plateaux and become free from temperature drift. 

A number of attempts have been made to obtain a reduction of the dead time. 
The discharge occurring in a counter operating in the Geiger region is stopped or 


quenched by positive ions produced during the discharge (Montgomery and 


Montgomery 1940). These ions form a sheath round the wire, reducing the 
field and preventing further multiple ionization. A second discharge is not 
possible until the sheath has moved to a critical distance from the wire, such that. 
the field is sufficiently restored. The time required for this movement of the ions 
has been identified by Stever (1942) as the dead time. 

Sherwin (see Smith 1948) has pointed out that, since the ion sheath is formed. 
within a few tenths of a millimetre of the wire, it should be possible to collect 


these ions on the wire, immediately after a discharge, by reversing the applied | 


field. Under these conditions, the times required for the positive ions to reach 
the wire would be relatively short and a considerable reduction in the dead time 
should be obtained. Following this suggestion Simpson (1944) and Hodson 


(1948) have used circuits which switched the wire potential a hundred or so 


volts negative with respect to the cylinder for a few microseconds after each 
count and then returned it to the working value. Both these workers used long 
counters up to 60 cm. Dead times of 20 to 30 psec. were obtained. 

If the mechanism by which the ion sheath spreads along the counter wire is 
considered, then another explanation of the reduction in dead times obtained by 
Simpson and Hodson is possible. Complete spreading of the ion sheath along 
the full length of the wire is a characteristic of the discharge of a counter operating 
in the Geiger region. Stever (1942) has shown that the ion sheath rarely spreads 
past a bead mounted on the wire in self-quenching counters. In a similar way, 
reducing the wire potential below the threshold of the Geiger region during the 
spread of the discharge will cause localization of the ion sheath. 

A consideration of the velocity of spread of the ion sheath (Wantuch 1947, 
Hill and Dunworth 1946) suggested that reversing circuits similar to those used 
by Simpson and Hodson, when used with long counters, should cause the wire 
potential to fali below the threshold of the Geiger region before the sheath has had 
time to travel far along the wire. The remainder of the counter would then be 
sensitive as soon as the wire potential is restored to its operating value. Thus the 
reduction of the dead time obtained with a switching circuit may be due to two 
causes: collection of positive ions or localization of the ion sheath. Some of the 
experiments described in §3 of this paper show that localization of the ion sheath 
does occur in long counters and that the explanation given by Simpson (1944) and 
Hodson (1948) of the reduction of the dead time obtained by them in their experi- 
ments was not complete. This conclusion (Collinge 1948) is in agreement with 
the work of Smith (1948), Hodson (see Smith 1948) and Elliot (1949). In this 
paper further experiments are described in which was determined the effect on 
the dead time of momentarily changing the wire potential. 


§2> THETCIR CWI 


The circuit which is shown in Figure 1 was intended for experimental purposes 
only. The cylinder of the counter is earthed and the wire is connected to the 
anode of the switching valve, V,,. The anode current of this valve is normally 
cut off by negative bias on the control grid. The cathode is returned to a preset 
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potential ’, which may be either positive or negative with respect to earth. Its 
anode is connected to a potential variable from 0 to 2,000 volts through a 10kQ 
resistance. Voltage pulses from the counter produced by ionizing particles are 
amplified, squared and fed on to the grid of V,, causing it to conduct. For the 
duration of these pulses the potential difference between the counter wire and 
cylinder is held at a value predetermined by V. It is convenient to call this 
potential difference V, the ‘switching voltage’. At the end of each pulse Vy, 
is cut off again and the wire potential rises to its normal value. 

The valves V, and V, amplify the pulses from the counter and feed them to 
the flip-flop formed by valves Vz and Vy. The ‘switching time’ is determined 
by the grid resistance of V, and may be varied from 1 to 10 psec. ; its amplitude, and 


Counter 


0-002uF 2-5kv. WKG. 
== 


Figure 1. Geiger counter switching circuit. 


hence the drive to V4, is determined by the cathode resistance Ol eVarur eo De 
circuit is ready to respond to a second pulse within 2 psec. of a previous pulse. 
This is an essential feature of the circuit and is ensured by the various D.C. restoring 
diodes. 

The value of the switching voltage was measured by connecting the counter 
wire to the Y-plates of a cathode-ray oscillograph. ‘The spot was normally 
biased off the screen by the high potential; each count caused it to be returned to 
the screen for the duration of the switching time and the switching voltage was 
read on a calibrated scale. This arrangement also enabled the drive to Vj, to 
be adjusted so as to obtain square pulses. 


§3. LOCALIZATION OF THE ION SHEATH 
(i) Preliminary Work 

The operation of a counter with the circuit described is as follows. If localiza- 
tion of the ion sheath occurs, the region of the counter, not involved in a particular 
discharge, will be sensitive to a second ionizing particle as soon as the wire potential 
is restored and the circuit is ready to receive a second pulse. Thus using an 
oscillograph, triggered by pulses from a counter operated with the circuit, it 
should be possible to observe counts occurring at all intervals greater than the 
insensitive time of the circuit. This insensitive time is the sum of the switching 
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time, the recovery time of the circuit and the time required for the wire potential 
to return to a point on the plateau where the counter pulses are large enough to | 
operate the circuit. a | 

To test this, several counters of lengths varying from 5 to 30 cm. and containing 
argon plus alcohol mixtures to a total pressure of about 11 cm. were used. Each | 
counter was connected in turn to the circuit and the anode potential of V,, adjusted | 
to the operating voltage for the counter. Pulses from the cathode of Viz ele used 
to trigger an oscilloscope, whilst the pulses occurring at the control grid at V, were | 
fed to the Y-plates. A gamma ray source gave a convenient counting rate. 
The switching time was varied from 1 to 10 psec. and the switching voltage from 
plus to minus 500 volts. 

Pulses were observed to occur on the cathode-ray oscillograph trace at all | 
time intervals 2 or 3 usec. greater than the switching time. The recovery time of 
the circuit was about 3 usec. ; this includes the time required for the wire potential 
to be restored. This experiment seems to confirm therefore that localization of | 
the sheath was taking place. 


(ii) Six-Cylinder Counter 
Direct proof that localization of the ion sheath was occurring in these long 
counters was obtained using the six-cylinder counter shown in Figure 2. The 
cathode consists of six, separately insulated, identical cylinders, contained in 
a ylass envelope. It was filled with 10 cm. of argon and 1 cm. of ethyl alcohol. 
The wire was connected to the anode of V,, and 1 MQ resistances were connected 
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Figure 2. Six-cylinder counter. 


Table 1 
Counts per minute at cathode No. 1 No. 2 No. 3 
Circuit inoperative ¥ 1124 1150 1131 
Circuit operative 918 97 43 


between each of the cylinders and ground. A discharge, occurring in a particular 
cathode, caused a positive pulse to appear across the associated resistance. This 
counter was arranged so that a collimated beam of gamma rays was directed, at 
right angles to the wire, on to No. 1 cathode. 

With a switching time of 1 psec. and zero switching voltage, the counting rates, 
at cathodes No. 1, 2 and 6, were determined by connecting an amplifier, dis- 
criminator and scale of 100 to each of the cathodes in turn. Readings were also 
obtained with the switching circuit made inoperative by disconnecting the 
condenser feeding the control grid of V,,. The results obtained are given in 
Table 1, which shows the average counting rate registered, at the three cathodes, 
after correcting for the background counting rate. 
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(iit) Two- Window Counter 

The counter shown in Figure 3 was used to obtain further confirmation in a 
similar way. It consists of a brass cathode with a mica window at each end, 
parallel to the wire. An RCA photo-multiplier tube was mounted opposite each 
window so that the light, produced during the discharge in that region of the 
counter, fell on the light sensitive surface. ‘The voltage pulses, occurring at the 
collecting electrode, were recorded on a scaler via an amplifier and discriminator. 

As before, the switching voltage was adjusted to zero and the switching time 
to 1 usec. Experiments were performed with either a collimated source of y-rays 
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Figure 3. ‘Two-window counter. 
Tables 2 
Counts per minute Window near source Window remote from source 
y-rays B-1ays B-rays y-rays B-rays B-rays 
Circuit inoperative 2051 3614 1687 1693 3736 1674 
Circuit operative 2087 3526 1568 220 42 16 
a SSS ae 
abies) 
—— EE ae 
Counter Six-cylinder Two-window 
‘Type of radiation y y y B B 
Distance between centres of detectors (cm.) 5 Des 62 32 32 
°% of discharges remote from source 10:5 4:5 1355 il 1 


directed at one end of the cylinder or with a beam of B-rays from **P passing into 
the counter at right angles to the wire, through an aluminium window. ‘This 
window, not shown in Figure 3, was in the brass wall opposite one of the mica 
windows. 

In Table 2 the number of counts per minute registered by the two photo-tubes 
are shown for both f and y-rays. The figures are average values and have been 
corrected for the background. 

(iv) Discussion 

The results of Tables 1 and 2 taken with the switching circuit operative are 
condensed in Table 3. The bottom line shows the percentage of discharges 
detected at various distances along the counters. These results show that 
localization of the ion sheath does occur. 

A higher percentage of the discharges was localized when f-rays were used ; 
this is explained by incomplete collimation and scattering of the y-rays. Even so, 
the results for the six-cylinder counter show that about 90° 4 of the discharges 


was confined to less than 5 cm. of the wire. 
BS2Z 
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§4. DETERMINATION OF THE DEAD TIME OF A LONG COUNTER 
(i) The Paired Source Method 


It was stated in § 3 that the dead time of a long counter obtained from atriggered 
oscillograph was equal to the insensitive time of the switching circuit. The dead 
time determined by this method is independent of the length of the counter 
provided it is long enough for localization of the ion sheath to occur. Thus with 
a switching time of 1 psec., dead times of 3 or 4 sec. were found for several 
counters of lengths greater than 5 cm. Although the result is consistent with | 
Korft’s definition of dead time, it is of no value for the correction of counting losses.. 
This is because the whole of the counter is not sensitive to further counts immedi- 
ately after the dead time. Regions of the counter where the ion sheath 1s localized 
remain insensitive for the normal dead time imposed by the ion sheath. The 
average losses depend on a number of factors which include the length of the 
counter. The true dead time must be obtained by a method giving the losses 
directly. 

A convenient way of determining the true dead time involves the use of two | 
sources. ‘The average counting rate is noted when the sources are presented to 
the counter, separately and together. 

If the registered counting rate is 7, then the true counting rate m, 1s given by — 


n 1-- : ou 
ee 1 
nm L+mr’ (1) 


where 7 is the dead time. 
From equation (1) Beers (1942) has derived a convenient formula for the 


dead time. Neglecting the background counting rate this may be written 


A Nya? 
2 2? 
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T= 


where A=n,+mn.—m,. and the counting rates ,, my) and n,, are-derived from 
N,, Nz and Nj, the total counts registered by the separate and combined sources. 
during the intervals ¢,, ¢, and t,) respectively. 

The error in the dead time can be calculated approximately. Suppose that 
errors in time measurements are neglected and the errors in the total counts. 
arranged to be of the same order by making N,=N,=N,,~WN say, then 
t, =t,=2t,,—t for roughly equal sources. 

Then from equation (2) it can be shown that 57, the probable error in 7, is. 
given by equation (3) where SN is the probable error in N. 


Oe aes 


The probable error P, due to statistical variations in the total number of counts. 
is obtained by putting IVN=2\/N 


, 2 t 
ser P= S| ee 
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from which it follows that 
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Then for chosen values of P, equation (4) gives the necessary counting period f 
for various values of t,5/N,». Some values of ¢ in minutes are given in Table 4, 
row 2, for P, = 3usec. 


Table 4 
ee eS Ea 
Counting rate 1. (sec~') 1,000 2,000 3,000 10,000 30,000 
Minimum counting period ¢ neces- 
sary to ensure P,<3 psec. (min.) 10 =< 2 202 <2 <2 
P,, the probable error in 7 fora=$% 
(usec.) 12 6 5 1 0-3 
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Errors due to uncertainties in the wire potential may often be neglected when 
finding the dead time of a counter operated in the normal way. ‘This is not true 
when a switching circuit is used which draws an average current which may be of 
the order of a milliampere, because of the difficulty of ensuring the same counter 
voltage at various counting rates. 

If the probable error in N, due to setting the counter voltage supply is 5N, 
then SN/N is the fractional error, « say. Then from equation (3) 


Pee lok eee (5) 


Row 3 of Table 4 shows the values of P, for « =4% for various values of Nijs/ ts. 


(ii) Experimental Details and Results 


The counter used consisted of a brass cathode, 32 cm. long and 2 cm. inside 
diameter, with a 0-2 mm. tungsten wire, and filled with 10 cm. of argon and 1 cm. 
of ethyl alcohol. In order to minimize the effect of positive ion collection the 
switching voltage was set to zero and the switching time at 1 psec. A scale of 
104, with a paralysis time of about 5 psec. was used to record the counts. ‘Two 
gamma-ray sources were used at a distance adjusted to give the desired counting 
rate. Readings were taken with the sources in several positions and the dead time 
determined as a function of the counting rate. A similar series of measurements 
was made with the circuit inoperative. The counting interval was always 
chosen from Table 4 so that the error due to statistical variation of the total 
counts could be neglected in comparison with «. 

The dead times obtained are shown in Figure 4 where they have been plotted 
against the counting rate. ‘The values of the counting rate, obtained with the 
two sources together, have been plotted along the x-axis after correcting for losses. 
Curve (a) was taken with the circuit not working and (5) with it working. 
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(iii) Discussion 

As the experiment was arranged so that the statistical errors could be neglected 
we need only consider errors due to setting of the H.T. The experimental 
arrangement enabled its value to be determined to about 5 volts. The slope of 
the counter characteristic at the wire potential used was 0:1°% pet volt, determined 
under the conditions for the dead time measurements. Thus in this experiment 
a=4%. Hence the values of P, applicable to the dead times obtained are not 
greater than those given in row 3 of Table 4. These errors apply to both graphs. 
The slope of 0-1% per volt is typical of the characteristic of a counter of this type, 
taken with a short paralysis time. 
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If the statistical dead time depends strongly on the counting rate, the value 
obtained from equation (2) loses its significance. Its meaning, in fact, is defined 
by this equation. This point has been discussed by Curran and Rae (1947). 


The values given in curve (a) Figure 4 are therefore somewhat artificial. The — 


results show however the large changes which occur when the counting rateis high 
and emphasize the difficulties of correcting for large losses. 
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Figure 4. Dead time as a function of counting rate for a long counter. 
(a) Circuit inoperative; (6) Circuit operative. 


The results shown in curve (6) obtained with the circuit operative show that 


not only is the dead time considerably reduced but its value is constant and 


independent of the counting rate within the experimental error. This is true up 
to counting rates of nearly 2 x 108 counts per minute. It would be reasonable to 
correct for losses greater than 10%, under these conditions. For a statistical dead 
time of 20 psec. the losses would be 10°% at a counting rate of 3 x 10° counts per 
minute. 

§5. EXPERIMENTS WITH SHORT COUNTERS 

In the experiment described in this section G.E.C. type G.M.2 counters were 
used; they are single ended argon-alcohol 8 counters. The tungsten wire is 
about 2.cm. long. Such a counter was connected to the circuit and a gamma-ray 
source used to provide a suitable counting rate. 

The dead time of such a counter was observed on a triggered oscillograph in 
the manner previously described. Dead times considerably greater than the 
switching time and dependent on the circuit parameters were obtained, as would 
be expected if the wire is too short for localization of the ion sheath to occur. 
To show that localization was not occurring, V, was used as a simple discriminator. 
The cathode was biased with a positive potential which was increased until 
triggering occurred only with pulses of maximum amplitude. This is a simple 
way of introducing a time delay to ensure that triggering occurs only when the 
ion sheath has had time to spread along the whole length of the wire. No change 
in the dead time was observed, showing that the ion sheath always spreads along 
the whole of the wire. 


{ 


i} 


Dead Times of Self-Quenching Counters 


As complete spreading occurs in these counters it was considered reasonable to 
bserve the dead time on a triggered oscillograph. The values obtained are of 
ualitative interest and are related to the values which would be obtained using 
he paired source method. 

The dead times obtained for one of the counters which had a dead time of 

105 psec. when the switching circuit was inoperative, have been plotted as a 
unction of the switching time on Figure 5. Several curves are shown for various 
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Figure 5. Variation of the dead time of a short counter with length of switching pulse. Vs is the 
switching potential. 


values of the switching time. For these results the anode load of V,, was changed 
from 10kQ to 50kQ. 
The curves show an expected initial decrease in dead time as the switching 
time increases. This is consistent with a progressive increase in the number of 
- positive ions collected on the wire. As the minimum dead time is not the same 
~ for all the curves and is very much greater than the switching time, it is concluded 
that part only of the ion sheath is collected. This is in agreement with the 
observed distribution of pulse amplitudes which shows an increase in amplitude 
from zero at the dead time to a maximum value. _ It is presumed that the decreased 
dead time corresponds to the time taken for the reduced ion sheath to move to a 
critical distance from the wire at which Geiger avalanches can occur. 

As the switching time is increased the rate of collection of -positive ions 
decreases. This is because it is more difficult to collect ions remote from the 
wire away from the very high field region. ‘The rate of collection decreases until 
the reduction in dead times obtained just equals the increase in the switching time. 
Thus the curves should have a minimum and a positive slope, approaching unity, 
for Jarge values of the switching times. This conclusion is in agreement with the 
curves shown. 

Some of the curves of Figure 5 do not include dead times for the full range of 
switching times. This is because for these conditions the cathode-ray oscillograph 


showed that very many spurious counts were occurring. 
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The reason why a reduction in dead times is obtained even when the switching 
voltage is negative has been discussed by Smith (1948). ‘The charge induced by 


thesheath on the wire leads to large fields near the wire sufficient to cause collection 


of positive ions. 
§6. COUNTER CHARACTERISTICS 

In general no significant changes in counter characteristics have been observam 
when using counters with the switching circuit. One counter had a somewhat 
longer plateau when the switching voltage was minus 150 volts. "Three commercial 
gamma-counters had plateaux with slopes which became progressively worse as 
the switching voltage was made more negative. All the counters tested produced 
larger numbers of spurious counts when the switching voltage was very negative. 

In order to reduce the chance of producing spurious counts, it is probably 


an advantage not to collect positive ions at the wire but to reduce the dead time 
of long counters by localizing the ion sheath. In this case a simpler circuit could 


be used to produce a very short negative pulse at the wire with an amplitude of 
300 or 400 volts. It would be desirable that such a circuit should trigger on very 


small pulses in order to ensure that the time available for the ion sheath to spread | 


is as small as possible. Not only is it desirable to limit the extent of the ion sheath, 


in order to obtain a short dead time, but also because the life of the counter is 


increased by reducing the number of positive ions involved in each discharge. 
Elhot (1949) describes some experiments in which the charge collected at the wire 
is shown to be reduced when a switching circuit is used and finds that under these 
conditions the life of the counter is prolonged. 
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ABSTRACT. When, in cold-cathode discharge tubes, electrode spacings and surface 
areas deviate from the ideal of infinite parallel planes, it is found that the breakdown voltage 
‘between electrodes increases. This is explained as a decrease in the effective » of the gas 
-due to loss of electrons and positive ions by diffusion to the walls of the container. The loss 
factor per unit potential difference is shown theoretically to be inversely proportional to the 
field strength and tube radius, although this is not fully verified by experiment. 

Paschen curves obtained experimentally for potassium and nickel cathodes in argon and 
jn a neon—argon mixture at low pressures are shown and, from these, values of y for the two 
‘cathode surfaces are obtained as a function of E/p). The apparatus used for measuring 
‘breakdown voltage is described. 

Earlier theories of statistical and formative time delays are extended to cover the case of a 
‘rising overvoltage and also the case where the primary electrons appear in bursts as with 
ionization by a-particles. The shape of the statistical distribution curve is an indication of 
whether the primary electrons have been produced singly or in bursts. ‘The overvoltage 
‘AV of breakdown due to formative lag and the rate of rise of the uniformly increasing 
applied voltage V aie found to bear the relation AV © \/(dV/dt). ‘This has been experi- 
mentally verified. 


§1. INTRODUCTION 


H1s subject has been treated by quite a number of authors and experimental 
work has been-carried out with various electrode macerials and gases. ‘The 
mechanism of breakdown is fairly well understood but the factor y which 

represents the secondary mechanisms has not yet been fully explored because of 
the difficulty of direct determination. In the following pages some possibilities 
care set down in this connection. 

When electrodes of finite dimensions are used in breakdown measurements, 
the special case of the similarity principle known as Paschen’s law does not hold 
as will be shown later and an approximate theoretical approach will be attempted. 


§2. INFINITE PARALLEL PLANE ELECTRODES 


The simplest breakdown gap from a theoretical point of view is that between 
infinite parallel planes giving uniform field strength. 

When all primary electrons are emitted from the cathode, the current density 
increases in the field direction according to the following relation (Penning and 
Druyvesteyn 1940, Rogowski 1940): 


[POS OS a” = (1) 
where V is the potential with respect to the cathode. 
If all primary electrons are due to ionization in the gas by an external agent, 
the current relation is 


fe rlexpeen Warley a) sole” 5 Rees ~- (2) 
In the steady state the positive ion current density at the cathode which 
corresponds to the electron current density of equation (1) is 


Jui =Jkel EXP {n(V, a Vo)} Sioa Mee hei an skates (3) 
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The probability that a positive ion which reaches the cathode causes emission 
of an electron may be represented by y,, and the probability that the electron will 
not return to the cathode by diffusion by y,,._ The probability that a positive ion 
will give rise to an active electron is thus: y;=Yiayip- 

Let there be @ photons created in the gas for each ionization and let a fraction 
fy fall on the cathode. The average probability that a photon will give rise to an_ 
active electron may be represented by y,,,, in analogy to the case of positive ions. 
The secondary emission factor due to photons is then (Penning and Druyvesteyn. 
1940) 


ua tyr — (Of (Vra¥ rb) 


The factor f, depends on the geometry of the gap and on photon absorption 
in the gas. In the present case, if absorption is negligible, f,=4. : 

Gis a function of the fundamental ratio E/p,) which is the main factor in discharge” 
at low pressure (Penning and Druyvesteyn 1940, Penning 1933). 


For metastable atoms an expression similar to that for photons will hold | 


(Penning and Druyvesteyn 1940): Fyn =(¢fn)(YmaYmb) is analogous to 6, 


both factors increasing with decreasing E/py. f,, depends on diffusion of the | 
metastable atoms and is therefore related to the geometry of the gap, the pressure: 
of the gas, and is also a function of time until the steady state is reached when _ 


metastable distribution between infinite parallel planes is (Rogowski 1940, 
Newton 1948) / 


(8) —ue Fey [exp (of Va— Vo} 1% —fexp (Ve Ve} =—DY sees (4) 
and 
c= Jae — Ae LE ate (5) 
eS al Vo Vo) exp in a= Vo) = 1 


The total secondary emission factor which relates the number of secondary 
electrons entering the discharge per second to the number of ionizations occurring 
in the gas per second is then (Penning and Druyvesteyn 1940) 


eam at teas zs ES trees (6) 


The factors y;., 7,. and y,,, are similar since they represent the eficiency of 
secondary emission. It has been shown by Penning (1930) and Oliphant (1929, 
1930) that the kinetic energies of ions and metastable atoms have little effect on 
these factors when the energies are low and of the order encountered in low pressure 
discharge. ‘The excitation energies are therefore responsible for secondary 
emission (Meili 1945) which means that the y, factors should be substantially 
constant over a large range of field strength and pressure. 

The factors yj, y, and yy, represent the effect of back-diffusion which has. 
been treated by several writers (Hertz 1927, Pose 1928-9, Thomson 1928, Loeb 
1939, Meili 1945). 

The general expression derived is (Pose 1928-9) 


ms oF sae (Vs “1 | 
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(This holds on the assumption that only elastic collisions occur in the distance: 
V/E from the cathode, which is approximately true.) It is evident from the 
expression that y, is a function of E/pp. 


| 
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It has thus been indicated that all the factors in equation (6) depend directly 

n E p, or are constants, except f,,, in the case of infinite parallel planes. If the 
contribution of the metastable atoms or the variation in f,, is not large, y will 
approximately be a function of Ep, alone. 

In the steady state f,, is a constant and if this state is approximately reached 
before final breakdown takes place y will be a function of E/po. 

The breakdown relation may be written in the form (Penning and 
Druyvesteyn 1940) 


1 1 
Ve=V,+—lo (145), A AN 8 
B 9 " g (8) 


The factors V and 7 have been determined experimentally (Penning and. 
Druyvesteyn 1940, Kruithof 1940) and are functions of E/p,) for pure gases. 

Thus the breakdown voltage V’, will be a function of E/po in the case of infinite: 
parallel plane electrodes. 

Since the field is uniform 


Pra 
eae oe ¢ Me ase 9 
AST © 


so that 7, will be a function of pod alone, which is a statement of Paschen’s law.. 
The minimum breakdown voltage corresponds roughly with the maximum. 
value of 7 (Penning and Druyvesteyn 1940) which for the rare gases is approxi- 
mately 1/3V;. 
The factor y varies enormously amongst various metals, being the greatest 
for the alkali metals and smallest for iron according to measurements made up to 


the present in the rare gases. 


§3, FINITE PLANE ELECTRODES 

When plane electrodes are used in measurements they have to be very large 
compared with their distance of separation to approximate to infinite planes. 
Loss of ionization products occur at the boundaries of the gap due to diffusion 
and the main cause of increase of breakdown voltage is the decrease in the effective 
7 as a result of the loss of electrons and ions. ‘The change in y will not be very 
great since only the factors F, and F), are reduced. When these factors become 
important at low values of E/p, this may not be true but changes in y do not affect 
the breakdown voltage so much, since it occurs in a logarithmic function. 

The greater part of the decrease in 7 is due to electron loss because of the large 
mean energy and consequent high rate of diffusion. 

In the breakdown or Townsend region at low pressure, space charge is 
unimportant (Penning and Druyvesteyn 1940) and it may be assumed that the 
diffusion of electrons and ions depends on their respective concentrations. 

At a distance x from the cathode the potential in a uniform field is V,= Ex 
and the potential difference over a distance dx is dV=Edx. 

The increment of electron current density due to ionization by collision is- 
njedV whilst the loss by diffusion may be represented by o,7 dV. 


Thus: dj,.=(n-oedV, 
feat EXP (9 — Fe Vaal oi, 1 Ses (10) 
The increment of ions is yj,dV and the loss oj,dV. Therefore 


dj; =—_ nJedV siz o,j,dV. 


Lf. G. fleymann 
Substituting for 7, and solving 
; é € f; 
jim Aexp [o(V—V)} ~ Us — exp {(-0,(V = Vo} 


Setting o=o, +0, and at the anode V,,=V,, j,=0, 


Ji jpe APT eM exp (qo) Vy Vo)} exp ((n- Ve Va}. 


The positive ion current at the cathode is the same as at V,,=Vg, 
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Breakdown will occur oe Vici =JIres 


UT [exp {(n-0)(V_—Vo)}-1]=1, 
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These are approximate expressions linking up the increase in breakdown 
voltage with the loss of electrons and ions due to diffusion. 

The factor o, is difficult to determine theoretically and the following derivation 
gives mainly a qualitative result: 

If the discharge is confined to the diameter of the electrodes of radius R by a 
glass envelope of slightly larger diameter and the electron concentration near the 
glass wall is n,,, then the number striking unit area on the glass wall per second is 


Gy =Nye,/4. ‘The average electron concentration in the tube is j,/eu, (7, is average 
-current density). 


Let n,, be related to this density by a factor a such that n, =aj,/eu,. ‘The loss 
per slenienel distance dx is then 


27Rqy dx =27R Tea. dx, 


Pag App E 
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“Therefore the loss factor per unit potential difference is 
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This is probably the main part of o and if o; is assumed to lead to a similar 
relation, the factor o/7 is approximately given by 
Veopoal 
= const. x on ER 

At constant E’/p, the only variables are E and R and with constant k, o/n should 
be inversely proportional to £ if the factor a does not vary. 

Mierdel and Steenbeck (1937) approach the problem differently by assuming” 
the average life of an electron or ion to be T= R?/(2-4)D. On this assumption 
o/ becomes inversely proportional to E*R?. 

In the experimental section the relative merits of the two derivations will be: 
discussed with reference to the actual results. 
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§4. MEASUREMENT OF BREAKDOWN VOLTAGE 


A diagram of the apparatus together with a sketch of the type of experimental 
tube used are shown in Figure 1. In measuring breakdown voltage, a steady 


400 v. 
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,_ to Discharge Diode 
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Figure 1. Apparatus for measuring breakdown potential and a typical experimental tube. 


applied through the potentiometer shown and the magnitude was indicated on 
voltmeter V,. On closing switch S,, the capacitance C charged through resistor R, 
and the potential across V, rose at the same rate, thus adding a further potential 
across the tube. When breakdown occurred, thyratron T fired, thus stopping the 
charging of C and applying an inverse voltage in series with V, and V, to reduce or 
extinguish the current through the tube, depending on the voltage. Vj, could read 
up to 240 volts and the inverse voltage was about 170 volts. 

The final value attained by V, together with V, is the maximum voltage across. 
the tube just before a discharge passes. We may thus define breakdown potential 
as measured by this equipment as the maximum voltage between the electrodes. 
prior to the negative slope transition region. 

The cathodes of the experimental diodes were illuminated sufficiently in order 
to reduce statistical spread and as will be seen later (Figure 8), the overvoltage at 
the rate of rise of V, of 40 volts per second is negligible. 

The experimental tube consisted of two mild steel electrodes enclosed by 


a close-fitting hard glass tube. One of the electrodes was made to slide on a ¢ in. 


} 
; 
potential about 100 to 200 volts less than the expected breakdown voltage was: 


= Lie We LLOVINGTTL 
diameter tungsten rod and could be moved by using an external magnet. The 
tube was evacuated and baked at 400° c. for some time and the electrodes were 
brought to a dull red heat by eddy current heating. In almost all the cases 
potassium-activation was used and a good, reproducible potassium layer could be 
deposited on the electrodes by vacuum distillation. 4 
Rare gases were used in the spectrally pure state and further purified by keeping 
in a vessel in which a continuous discharge was maintained between potassium- 
activated electrodes. Condensable vapours were removed in a double liquid 
-oxygen trap andas a last safeguard a potassium-activated diode passing a continuous 
‘discharge was permanently connected in the inlet to the main tube. 
Three sizes of tube (dimensions in centimetres) were used: 


Electrode Tube Maximum electrode 
Tube : ; ; 
diameter diameter separation 
1 1°5 1:6 Ds 
2) 2°8 3:0 2-6 
3 6-0 6:5 8°75 


Breakdown voltage was measured as a function of electrode separation at 
various pressures in argon and in a mixture of neon and argon which contained | 
6°6% argon (potassium-activated cathode). 

Some representative curves for tube 2 are shown in Figure 2. Breakdown 
voltages are plotted as functions of pod and it is evident that the curves at the lower 
pressures do not follow Paschen’s law. 

In Figure 3 the curves of breakdown voltage in argon are plotted for all three 
tubes as functions of E/p), assuming- uniform field distribution. As was to be 
expected, the differences between curves are greatest for the smallest tube and 
least for the largest, for the same differences in pressure. At higher pressures the 
‘curves tend to merge in all cases and the curves for the various tubes show good 
‘correspondence at the pressures at which the losses by diffusion become small. 

In order to determine the effect of the nearness of the glass wall to the edges of 
the electrodes, a tube with 1-5 cm. diameter electrodes and a 3cm. diameter glass 
envelope was constructed. The values of breakdown voltage were not greatly 
different from those in tube 1 at the lower pressures, and it was therefore decided 
that the glass wall had only a small effect. 

In order to check the validity of the theoretical expression for o/7 given in 
equation (15), the curves in Figure 4(a) were drawn, using equation (13) to deter- 
mine o/n. ‘The theory predicts that o/7 should be inversely proportional to £ 
at constant E’/py and it has been found that there is reasonable correspondence as 
shown by the experimental points in Figure 4(a). 

These points were calculated at various values of E/p) (from 50 to 
200 volt/cm.mm.), and it seems that for each tube a single factor of porportionality 
will represent all the points satisfactorily. The correspondence is not exact and 
partly due to experimental and graphical errors. Furthermore the theory is only 
‘approximate. 

The correspondence amongst the three tubes is not so good, the theory 
predicting that the slopes of the lines in Figure 4 (a) should be inversely proportional 
to the electrode diameters, which is not the case. It is probable that there is a 
‘further factor which has been overlooked in the theoretical derivation or it may be 
that the factor a in equation (15) is a function of R. 
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Figure 2. Breakdown potential curves with Figure 3. {Breakdown potential curves as 
electrodes of 2:8 cm. diameter at various functions of E/p) in pure argon. Electrode 
pressures : diameter : 

' {a) Pure argon; (b) 93-4% neon, 6:6% argon. (a) 1:5cm.; (6) 2:8cm.; (c)" 6-0 cm. 
: 
: 
: 


S\5 


0 0:01 0-02 0-03 0:04 
\/E (Z in v/cm) 


Figure 4. Experimental points representing o/7=AV/(AV + V). 


(a) Argon: potassium-activated electrodes (electrode diameters in cm.). 
(b) Argon: nickel electrodes. 
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There seems to be a further objection to the theoretical expression in that the 
factor o/7 for nickel electrodes of 6cm. diameter derived from Figure 4(d) does | 
not correspond to that for potassium-activated electrodes. 

The theoretical approach according to Mierdel and Steenbeck (1937) appa | 
even less than the theory derived in this paper and cannot therefore explain che | 
increase in breakdown voltage. 

Figures quoted by McCallum and Klatzow (1934) for increase in breakdowil | 
voltage between nickel electrodes in argon are much higher than those found in | 
this series of experiments and must be due to an additional cause not present here. | 

The increase of breakdown potential with small electrodes and at low pressures. 
makes it necessary to use large electrode diameter, small spacing and high pressure 
to obtain a true Paschen curve. The values of breakdown potential measured. 
on the tube of 6cm. diameter will correspond fairly well with the values between 
infinite parallel planes except perhaps at high values of pod in argon where the 
pressure was about 50mm. Hg and the separation up to 3°75 cm. 

The measured curves for nickel and potassium over an extended range of pod 
are shown in Figure 5. The curve for potassium is easily reproducible and on 
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Figure 5. Experimental breakdown potential curves in argon with potassium and nickel cathodes. 
Inset : y-curves as functions of E/p 9, for potassium and two different nickel cathodes. 


using each electrode as cathode in turn, very close correspondence between the 
breakdown potentials occurred. 

The nickel was in the form of an electro-deposit on the steel electrodes and 
although the treatment of both electrodes was exactly similar, the breakdown 
potentials with each as cathode in turn differed appreciably, 
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The factor y was calculated from the breakdown relation in equation (8), 
though high accuracy is not attained by this method. 

The values of Vy and 7 were taken from Penning and Druyvesteyn (1940) and 
ruithof (1940). Figure 5 (inset) shows y-curves for potassium and nickel as 
unctions of E/p). Curve a for nickel corresponds with the pod breakdown voltage 
urve shown. Curve bd is for the other electrode and merges with a at higher 
alues of Epp. 

The curve of y for potassium shows a form similar to that to be expected if the 
, factors were predominant, the other factors being approximately constant. At 
ery low values of E/p, the curve tends to rise again as has been found for other 

aterials by others (Penning and Druyvesteyn 1940) and is probably due to the 
increasing importance of @ and ys as E/p, decreases. 
The curves for nickel show that the metastable atom and photon action become 
important below a value of £/p, of 100 volt/cm.mm, and to adifferent extent for the 
two electrodes, depending on their sensitivity to this action. Again this confirms 
previous results (Penning and Druyvesteyn 1940), but the full explanation is still 
bscure. 

This illustrates the fact of the difficulty of the interpretation of y and the 
importance of the purity of the cathode material. 


§5. TIME DELAY EFFECTS IN BREAKDOWN 


There are two main causes of time delays in breakdown which will be treated 
separately and then combined. The delays involved when constant overvoltage 
is applied have already been dealt with by others (von Laue 1925, Zuber 1925, 
Schade 1936-7, Hertz 1937) and will serve as a basis for relations to be derived for 
‘the case of linearly increasing overvoltage. 

The statistical time lag is due to insufficient primary electrons and formative 
lag comes about as a result of the time taken for positive ions to move to the cathode 


from their points of origin. 


, §6. STATISTICAL TIME DELAY 

_ Primary electrons, which are necessary to initiate breakdown, appear in a 
random fashion due to some external agent and therefore introduce time delays, 
even when the applied potential is adequate for breakdown. Let formative time 
delay be neglected for the time being. Von Laue (1925) has shown that the 
probability that breakdown will take place after a time ¢, measured from the instant 
of applying an adequate potential, is given by 


a (=exp(~4 le Wat). ror (16) 


In the case of plane parallel electrodes, Hertz (1937) has shown that W, the 
probability that one electron will initiate discharge, is related to the factor M as 


follows: 
log 1-W 
M=yfexp (1(V.—Va)} = EE. 


At the normal value of breakdown potential M=1 and W=0. 
With constant overvoltage, i.e. when M is constant and greater than unity, 
W is constant and the probability of breakdown becomes 


filé)=exp (—G_.Wt,) =exp(—te/to). ete (17) 
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This relation has been verified experimentally by Zuber (1925). 

In order to extend the theory to the case of rising overvoltage, it 1s necessary 
to make some approximations : 4q 

If y and W are small at small overvoltage AV, then 


1 s 
M=- pics (l— W)=1+3W. 


+ 


Using a method of Schade (1938), it follows approximately that 
M=y [exp {(ap +Ax)(d—dy)} — 1] | 


=y[exp {ap(¢—do)} —1] +y exp tan(d— d)} Se ABd dy) 


lo 
=1+ a, AV. ss ap eeeee (18) | 
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Therefore the probability of breakdown after a time ¢, with the applied potential 
rising linearly from the normal breakdown value, is 


So(t) = exp {— 2ge(OW/ot)t.?} =exp {—2(ts/tm)?}- sven (20) 
The probability of breakdown occurring in an interval between ¢, and ¢, + dt, 
is given by 


fi dis = exp {iMG la) Alt), ae (21) 


with a maximum occurring at ¢t,=t,. 

The functions f, and f, are represented in Figure 6 in terms of their characteristic 
time constants, ft) and ¢,, respectively. 

Function f, is only valid for small values of overvoltage and cannot be applied 
where the overvoltage becomes large. It is further restricted if the electrons do 
not occur singly but in 5 bursts of r electrons each. The factor g, must then be 
replaced by br and Wby W,=1—(1—W)’. Whilst W (and therefore AV) is very 
small this does not differ much from f, since b(@W,,/dt) =br(dW/ot) =4,(dW/dt), 
but W,, soon approaches unity and the function goes over into a form similar to f,. 
A better approximation at larger overvoltages may be made by assuming W to have 
the form W=1-—exp(-—t/T), where T is approximately determined by finding the 
value of AV, which corresponds to W=1—e-}, giving T=AV ,/(dV /dt). 

The probability of breakdown for electrons occurring singly then becomes 


faa(t) =exp [ ae qetts = iG — €xp ( a> z./ T))}]. HORS OG (22) 


If the average number of electrons formed simultaneously is r with b bursts per 
second, 


fart) =exp[—b{t,—(T/r)\(1—exp(—rt/T))}]. (23) 
The probabilities of breakdown in any particular interval may be found by 
differentiation and it is found that maxima occur at the following values corre- 
sponding to equations (22) and (23) respectively : 
Agel: 
> t,/T=log ——— ++ 
Fan ctaeesh 85, Tres 4ecly. see (24) 
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Some general curves representing Soa and f;, are shown in Figure 7(a). In 
hese ¢,’=brt=2 is constant and r is varied. For any other value of q,T the 
robability is just raised to a proportional power. 

It is evident that already with bursts of 10 electrons at a time the curve fy, 
pproaches closely the exponential curve exp (—bz¢,). 

Thus the transition is fairly rapid from one shape of probability curve to the 
ther. The form of the statistical distribution curve will indicate the way in 
hich the primary electrons are produced. Where these electrons are due to 
mission from the cathode, they are independent of each other and therefore the 
distribution represented by fs is to be expected. 

When the primary electrons are produced by «-particles, they will occur in 
roups, each particle giving rise to a whole train of electrons which are simultaneous 
for most practical purposes, so that /;, should be applicable. 
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Figure 6. Statistical time delay probability curves /, and f, and curves of combined statistical and 
; formative time delays. 


Ionization by other types of radiation does not usually give rise to a closely 
spaced group of electrons and the distribution will tend to approximate to the case 
of single electrons. 

Electrons which are produced in the body of the gas are less effective than 
electrons emitted from the cathode. The probability function W is based on the 
latter and in the case of gas ionization it may be shown that W must be divided by 
n(V —V,) to compensate for this fact. 

Some interesting experiments have been carried out in order to determine 
statistical delay time when a discharge tube was kept in total darkness away from 
radioactive sources. ‘The oniy external source of ionizing radiation could have 
been cosmic rays and the very feeble natural radioactivity of materials. 


Cc-2 


36 F. G. Heymann 


The experimental tubes used were small sealed-off diodes with potassium 
activated cathodes and argon filling at 18mm. Hg pressure. These were of the | 
parallel plane electrode type, electrodes spaced 0-06 in. apart and }in. in diameter, 

It was found that primary electrons far in excess of those to be expected from 
the above-mentioned causes were present in the tube when the breakdown | 
potential was measured at intervals of three seconds. On increasing the interval | 
between measurements to 90 seconds, the primary electrons were much less, _ 
showing that the short discharge occurring at each measurement had some 
influence on the subsequent time delay. Between measurements the inter- | 
electrode field was reduced to zero. 

Representative curves are shown in Figure 7(4). The distribution curve for 
measurements at 3 second intervals yields a calculated rate of primary electron 
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Figure 7. 


(a) Statistical time delay probability curves /3, and fs). 


(6) Experimental statistical distribution curves with different intervals 
between measurements. (Rate of voltage rise—40 v/sec.) 


production of about 5 per second. With 90 second intervals the rate is about one - 
electron in two seconds, whereas the number which will be produced naturally is - 
of the order of one per minute. 

This effect may be attributed to continued emission from the cathode after © 
extinction of a discharge, as Paetow (1939) has shown. He found that this. 
continued emission was still appreciable after 50 seconds and that it was probably - 
due to small specks of impurity on the cathode which acquired charges during: 
discharge and afterwards caused electron emission due to the high fields set up. 
It is probable that even on high purity potassium such specks may exist and give) 
rise to this persistent emission which was still noticeable when the interval between | 
small measurement discharges was 90 seconds. | 

Paetow found that the emission was very large immediately after stopping a. 
discharge, partly due to metastable action which became negligible after about | 
a second and thereafter due to this after-effect only. | 
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§7. FORMATIVE TIME DELAY 
When statistical time delay is reduced to negligible proportions by sufficient 
primary electrons, the discharge will not take place instantaneously on application 
of sufficient overvoltage, since it takes a definite time for the current to increase to 
the point at which conduction is appreciable. 
This formative time delay has been dealt with by Schade (1936-37) for the 
case of constant overvoltage, yielding the following expression : 


ig=t, ca log (RO pe. ee (26) 
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The main cause of the delay in current increase is the average transit time of 
positive ions from their points of creation to the cathode. The primary electron 
current density at the cathode is j,) and the final cathode electron current density 
jxe and their ratio which is included in the logarithmic term has a minor effect on 
the formative time delay. 

The expression for formative delay has not previously been derived for the case 
of rising applied potential and will be developed below. 

By the method used before in equation (18) it is possible to show that 

Oo da dV 
M—1=57AVy= a6 dp 

For small overvoltages (M — 1) will be small and the time delay will be approxi- 
mately inversely proportional to the average value of (M—1), according to 
equation (26). 

: If the overvoltage rises linearly, the average overvoltage will be half of the final 
value reached 
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and from equation (26) 
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The terms in the logarithmic function can cause very little variation so that it 
follows approximately that 


tepoc(dVJdt)-#,  AV,,co(dV [dt saenee (29) 
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On deriving the expression more rigorously from Schade’s (1936-37) ditfer- 
ential equation, the result is 


1 Uke oud NANT tN 
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There is not much difference when compared with equation (27) which may 
therefore be used. 
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The limiting formative time delay when (IM —1) tends to zero becomes 
(Schade 1936-37) 


LEO —Cilc ice a eee eee (31) 4 
M->1 


The average distance through which positive ions move to the cathode may be 
found by integrating over the whole discharge gap, yielding 


e*4 1 
ae =a a al Ci CO A RL ed mp nye cere (32) 


This is for plane parallel electrodes so that the average transit time is 
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In the case of a concentric diode in which the anode radius R, is much smaller 
than the cathode radius R,, the average transit time is approximately 
fe R, ie R, 
' 2K;E,/po’ 


where E, is the field strength at the cathode. 

There is some experimental confirmation of the equation (28) in that the linear 
portions of measured curves have slopes which correspond well with the calculated 
values. 


20 


200 


Figure 8. Experimental formative overvoltage curves at various values of po and d. ‘The linear 
portions correspond well with the theoretical calculations. dV /dt expressed in v/sec. 


Figure 8 shows typical curves from which it is evident that in general the initial 
portions of the curves do not conform to the theory but that the curves have a 
constant slope over most of the range of rate of voltage rise. 

The expected slope of an experimental curve is AV,,/s/(aV /dt), from 
equation (28). 
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Assuming log {1 +(M—1)ayJie/Jxo} to be approximately constant at 18 and 
aking 0x/0E from Figure 9 (Kruithof 1940) the following table may be drawn up for 
argon and 6cm. diameter potassium-activated electrodes to compare experiment 
with theory: 


Ne-A —E/py (V/cm.mm.) 
60 


0 20 40 80 Ke) 
0-04 
Ye 
oF 
0-02 
0 
0 100 200 300 400 500 
Argon—E/p) (¥/cm.mm) 
1:0 T 
N 
No 


0 0-1 0-2 0:3 0:4 0:5 
¢ (msec) 


Figure 9. (a) The factor @a/0F for argon and for a mixture of 93:4% neon with 6:6% argon. 
(b) Experimental distribution curve showing the effect of combined statistical and 
formative delays. 


Experimental Theoretical 
Po d (5 Slope Slope 
4°56 0-104 Ori 0:023 0-017 
0-195 0-53 0-031 0-028 
0-405 1:74 0-050 0-049 
0:8 4-95 0-086 0-083 
9:5 0-1 0:29 0-025 0-020 
0-2 0:89 0-035 0-035 
0-41 2:74 0-065 0-064 
8:65 0-1 0:27 0-025 0-020 
0:2 0-84 0-036 0-034 
0-41 2°58 0-050 0-061 


0-81 7:05 0-114 0-107 


Mobility of argon positive ions at 1mm. Hg = 2:1 x 108 cm/sec/volt/cm. 

The theory thus finds some justification in practice except at small overvoltage 
where more factors seem to be present than those which have been taken into 
account. The order of overvoltage calculated is correct and the theory may 
therefore be applied with reasonable accuracy. 

This formative delay is not a variable quantity in a fully symmetrical discharge 
gap and may therefore be allowed for in cases where it becomes important. 
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§8. COMBINED STATISTICAL AND FORMATIVE DELAYS 


2 


If average statistical and formative delays are of the same order, an expression ! 


may be derived in order to take this into account for rising potential. 


With statistical time delay, a certain time elapses before the discharge is 
initiated and therefore the overvoltage will reach a value AV, before the formative : 


period begins. 
Let the increase of overvoltage during the formative period be AV,. Then. 
the average overvoltage during this period is AV, + 4AV,. 
The formative time delay is therefore given by equation (26) as 


1 Zh; tae 
ae AV,+4AVy 0a /0E log ul (fl L)avJielIxos ’ 
av, ___AVee 


Or a BY AAV 


where Vg is the formative overvoltage when V,=0. 
The solution to the above equation is 


AVa= —AV, +{AV2+AVa 9 eee (35) 
The total overvoltage is | 
AV=AV,+AVy={AV 2 +R, ts ae (36) 
f= [Ee tig t=O 2 eee (37) 


The statistical probability curves f, and f, are modified as shown in Figure 6 
when the formative time delays are of the same order as the statistical delay time 
constants ¢) and ¢,, (see equations (17) and (20)). 


A typical statistical distribution of points for combined delay times is shown in — 
Figure 9(b), as measured on a small concentric diode containing 15mm. argon | 


and about 1 wgm. y-equivalent radioactive material to provide primary ionization. 
The measured formative time delay at a rate of voltage rise of 3 x 10! Ko 


was about 0-25 millisecond. The solid curve was fitted at the ordinate e? to — 
correspond to a modified f, probability curve with ty) =2¢,,.. The broken curve | 
was fitted at the ordinate e-' and corresponds to a modified f, probability curve with | 
tg =to. No better fit is possible with this latter type of distribution and it thus | 


seems as if the /, type of distribution given in equation (20) gives a very good fit. 


This is however not what would be expected if the «-particles created groups of 


electrons containing about 50 electrons each as they should do in this case. 

The number of primary electrons occurring per second on calculating with the 
J, type of distribution is about 4 x 104 per second, which is of the order of the 
number of «-particles to be expected per second. 

Theoretically the curve should follow the f, type of distribution which in this 
case does not seem to be confirmed. Other tubes show similar curves which 
means that there may be a further factor to take into account. 
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: The Behaviour of Multiple Circuit Magnetrons in the 
Neighbourhood of the Critical Anode Voltage 


By W. E. WILLSHAW anp R. G. ROBERTSHAW 
Research Laboratories of The General Electric Company, Limited, Wembley, Middx. 


MS. received 30th Fune 1949 


ABSTRACT. It is pointed out that the mechanism of operation of the multiple circuit 
magnetron oscillator in the region of minimum magnetic field and voltage, where the 
efficiency is commonly assumed to approach zero, should approach that of an oscillator of 
the travelling wave tube type, providing that a cathode of suitable size is used. Useful 
efficiencies should thus be obtainable under these conditions. 
Details of experiments are given in which an electronic efficiency of 12°% was obtained 
at a wavelength of 3 cm. at values of magnetic field and voltage several times lower than those 
used for high efficiency operation. The mode of operation was determined by the value of the 
magnetic field, a given mode being maintained over a range of magnetic field of the order of 
8%. The anode voltage was about 70%, of the critical value. 
The experimental results generally support the hypothesis and suggest that the minimum 
voltage regime should be of extreme importance for work at the highest radio frequencies. 


SiN CDR OD WiC TNOIN: 
HE main characteristics of performance of the high efficiency multiple 
circuit magnetron oscillator are well known. Operation in one mode 
starts at mimimum values of magnetic field strength and anode voltage, and 
continues with increasing efficiency as these parameters are increased. In the 
higher efficiency region it is assumed that electrons interacting favourably with 
the oscillating field leave the main electron cloud, which extends only part of the 
way from cathode to anode, and being gradually focused into tight bunches 
finally arrive at the anode with a velocity which is the same as that of a rotating 
wave component of the anode field. On this assumption a very good approxt- 
mation to the observed efficiency is derived directly, i.e. y=1—}mQ?a*/eV- 
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where is the rotating wave angular velocity, a the anode radius, and V, thei 
anode voltage (Willshaw and Rushforth 1946). Change of mode, giving a change) 
in Q may occur at any value of magnetic field strength as V,, is changed in order} 
to change current. For the lower efficiency region, where the anode voltage any 
magnetic field strength are much lower, the main electron cloud extends muchl 
nearer to the anode, so that electrons leaving it do not have time to be formed 
into bunches before reaching the anode, and to acquire the velocity of the rotatings 
wave component. ‘Thus the observed efficiency is much lower than that given) 


by the above expression. 


§2, BEHAVIOUR NEAR THE CRITICAL VOLTAGE 


At the limiting lower values of magnetic field and anode voltage the electron 
cloud extends quite close to the anode surface, and there is experimenial evidences 
to show that in the absence of oscillations electrons may describe circular paths 
which are coaxial with the anode cylinder (Wassermann 1948). Now the electro 
velocity at radius 7 is eH/2m(1—b6?/r?), H being the magnetic field strength,| 
and 6 the cathode radius. If the radius of the cathode is not too large relative to; 
that of the anode, this velocity will vary very little in the neighbourhood of the¢ 
anode. The variation in electron density in the neighbourhood of the anode iss 
similarly small, being proportional to 1+(b/r)*.. Conditions are thus similar to¢ 
thcse in a travelling wave tube (Kompfner 1947), in that there is a stream of 
electrons of given velocity and density passing across the mouths of severali 
cavities and therefore capable of interacting with the fields due to oscillations 
in these cavities. By direct analogy there is then the possibility of oscillations$ 
being generated when the electron velocity close to the anode is equal to theé 
velocity of the wave generated in the anode circuit. This condition will obtain: 
when (eH/2m)(1 —}?/r,,?) =Q=(2eV,,/m)! where r,, is the maximum radius/ 
reached by the space charge. A change in mode of oscillation corresponding to 
a change in Q should result when H or V,) are correspondingly changed, this: 
being in distinct contrast to the behaviour in the high efficiency region. Wed 
shall call this the ‘minimum voltage’ regime of operation. 


§3. EXPERIMENTAL PERFORMANCE IN THE MINIMUM 
VOLTAGE REGIME 

Some experiments have recently been carried out in the minimum voltages 
regime with a magnetron having 28 cavities of the ‘rising sun’ form (Hollenberg, * 
Kroll and Millman 1948) with a circuit ratio of 2-:25/1, anode diameter of 16 mm,, . 
cathode diameter of 3 mm., and wavelength of 7 mode of 3-31 cm. Figure 1. This: 
valve, fitted with a much larger cathode (115mm. dia.) had previously been: 
designed and used for high power operation with an anode voltage of the order: 
of 15 kv., magnetic field of 2,500 oersteds and efficiency of 20°. Figure 2 shows: 
the relationship between anode voltage and mean current for the minimum 
voltage regime for a range of magnetic field values, and contours of constant! 
power output are included. The valve was operated with a 50 c/s. voltage applied | 
to the anode through a 5,000 ohm series resistance so that current passed for 
approximately one fifth of the time. The corresponding wavelength of operation 
is shown together with the mode number n of the wave component with which: 
the electron is interacting. The mode number is physically the number of 
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Figure 1. Rising sun magnetron for 3:3 cm. wavelength, and water load. 


Anode Voltage \ KV.) 


Mean Anode Current (ma) 


Figure 2. Variations of anode current and mode of operation with anode voltage and magnetic 
field strength. (H—=magnetic field strength in oersteds, A=wavelength in cm., ny =funda- 
mental mode number, n=wave component mode number.) 
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times a wave is repeated round the circumference of the anode (Willshaw e¢ al 
1946). The mode number z, shown in the Figure is the fundamental mode of 
which u is a component. It can be seen that the mode of oscillation depe dd 
-entirely in the magnetic field, and oscillations are maintained in a given m dd 
over a range of magnetic field of the order of 8°, a maximum efficiency of 6%) 
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Figure 3. Relation between the observed and calculated anode voltage and magnetic field strength. 
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Figure 4. Relation between voltage, efficiency and magnetic field strength for minimum voltake 
and high efficiency operation. 
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ing obtained in the mode n=14. The electronic efficiency in this mode, 
timated from the ratio between power dissipated in the load and in the valve 
reuit, is 12°. In contrast to this, when operating as a higher efficiency pulsed 
ve, at an anode voltage of about 15kv., oscillations are obtained in the mode 
=14 (the 7 mode) over a wide range of magnetic field strength, as is usual for 
ves operated in this regime. Figure 3 shows the observed operating voltages 
otted as a function of magnetic field, together with the values calculated for the 
ifferent modes of operation on the hypothesis of the minimum voltage regime, 
suming that the outside of the electron cloud reaches the surface of the anode 
that r,,=a. These calculated voltages are plotted as a function of the magnetic 
elds at which oscillations are observed and are rather lower than the observed 
oltages. The critical voltages are also indicated and it is seen that operation 
es place at voltages which are roughly 70% of these. Finally Figure 4 shows 
jagrammatically the relationship between the various parameters for the high 
ficiency and minimum voltage regimes of operation. The dotted extensions 
f the curves of high efficiency represent the values calculated from the simple 
xpression mentioned earlier, the practical values lying much below these and 
alling to zero at a magnetic field greater than that at which ‘minimum voltage’ 
oo occurs. 
§4. FIELD OF APPLICATION 
It is evident from the results which have so far been achieved that useful 
“ficiency may be obtained with a multiple circuit magnetron operating like a 
-ravelling wave tube in the minimum voltage regime. Values of anode voltage 
and magnetic field are several times smaller than for the high efficiency regime. 
The size of the valve for a given voltage is thus a maximum, and this fact and the 
low magnetic field involved make this mode of operation of extreme importance: 
for work at the shortest wavelengths. 
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Measurements of the Reflection Coefficient of he 
Water at a Wavelength of 8:7 mm. 


By Di 4G. KIELY: 


Royal Naval Scientific Service 


Communicated by L. H. Bainbridge-Bell ; MS. received 25th April 1949, 
and in amended form 3rd August 1949 


ABSTRACT. An account is given of measurements (carried out in March 1947) of the 
reflection coefficient of water at a wavelength of 8-7 mm. over a range of angles of incidence, 

The method employed is to measure the relative field strengths of the direct waves and 
waves reflected from a trough of water using free-space propagation and high-gain aerials. 

The following electrical constants of water have been computed from the measured 
results (water temperature 11:1°c.): refractive index 4:40-+0-24, dielectric constant 
10°86+2-21, absorption coefficient 2:91-+0-06, conductivity/frequency 12:82+0-42, 
Brewster angle 79°. | 


ol ENR ODiCiCa TON; 


HIS note describes some experimental work carried out at a wavelength of | 

8-7 mm. to determine the reflection coefficient, conductivity, dielectric : 

constant and absorption coefficient of water. The reflection coefficient - 
was measured for different angles of incidence, and from these results the other : 
constants were calculated. At this wavelength the problem of measuring the » 
reflection coefficient was considered more optical than radio in nature and the 
method employed was to measure the relative field strengths of direct and reflected 
free-space waves using high-gain (narrow beam width) aerials and a trough of 
water. 

§2. DESCRIPTION OF APPARATUS 


The experimental apparatus consisted of a transmitter and receiver, mounted 
on metal racks, which allowed the heights of the aerials to be varied, above a 
trough of water placed between them. The horizontal distance between the 
receiver and the transmitter was maintained constant and the total variation in 
height of either aerial above the water surface was approximately 4 feet. The 
trough dimensions were 8 ft. 9 in. by 3 ft. The aerials and radio-frequency 
circuits were mounted on tilting platforms attached to the racks: this allowed 
them to be pointed down at the water surface or directly across it. The beam 
width was sufficiently small for the direct and reflected rays to be received 
independently. Polarization was vertical. 


§3. EXPERIMENTAL PROCEDURE 


The horizontal distance between the receiving and the transmitting aerials was 
kept fixed. 'The heights of the two aerials were adjusted so that for each angle of 
incidence the aerials were directed at the centre of the trough. The experimental 
procedure was extremely simple and straightforward. The power was switched 
on, 30 minutes was allowed for the klystrons to become warmandstablein frequency, 
the transmitter and receiver were tuned, and the ceceiver gain was turned to a 
suitable level. For each angle of incidence a direct transmission across the trough 
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s received, and the receiver output, V,, and the transmitter power monitor 
ding noted. The aerials were then set to the required angle of incidence, the 

ected wave received, receiver output V, noted and the transmitter power 
onitor reading checked. This was repeated several times for each angle of 
cidence and the average value of V, and V, obtained. Only very small variations 
V, and V, were noted in these repetitions. A correction was made in the value 
V, for the difference in path length for direct and reflected waves using the inverse- 
stance law to give Vy; the ratio V/V, was then equal to | R’ +7R"|, the modulus 
the reflection coefficient, R’ and iR" being the real and imaginary components. 
he water temperature during the measurements was 11-1°c. The experimental 
sults are shown in the Figure as a function of angle of incidence. The point 
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e Experimental Points 
A =8-7mm. 
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t 87-5° was obtained at a wavelength of 6-35 mm. by reflections from the sea 
uring another experiment. This point has not been used in the following 
heoretical work. 

An experiment was carried out to check the validity of the inverse-distance law 
of propagation for the distances involved in the above work. This experiment 
showed that the inverse-distance law was valid for the six measurements from 
37°5° to 76°, but was not strictly valid for the two measurements at 19° and 30° 
nejidence. These two observations have not been used in the following 


calculations. 


§4. THEORETICAL ANALYSIS OF RESULTS 


The symbols used are as follows : »2=complex dielectric constant of water, 
n=real part of refractive index of water, «=real part of dielectric constant of 
water, « =absorption coefficient of water, f=frequency of radiation, R=R’ +7R" 
=complex coefficient of reflection of water, 9=angle of incidence, ¢ =complex 
angle of refraction, o =conductivity of water. 

It is suspected that «>1 and o/f>1. Now 
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If sind=-1sin 8, the reflection coefficient is 


_ | ecos@—cos¢d 2 = | 
[Rss weosi feos | ee ( a) | 
Now, as there is @ priori reason to suspect that |u?| =|e—7(20/f)| >1, an 
since sin@<1, 
sin? 6 sin? @ | 
cosp=/(1- 7D )=1-3 2 ee Seen fle i> (2 5) 


and the error in adopting this approximation will be at the most of the order of 2% 
Hence 
| u2 | cos?6 + 1 —2ncosé 


2 = ea ep 
Et | 42 | cos?64+1+2ncos6° 
eR |e aye? heos20 et 
Let wee SR? | Opcus Gn 22 eee 5 (3) 
ae: ys On+ binge 2 oT eee (3a) 


where x=cos0; a=|,.?|/2n; b=1/2n. In order to determine the best values of} 
a and 6, the method of least squares was adopted with suitable weighting factors. 
This gives the values a=3-163 + 0-065; 6=0-114 + 0-006. 

Having found a and 4, the physical constants for water were obtained as: 
follows: 4 


n=1/2b = 4-40 +0-24, 
€= 1/26? —a/b =10-8642.21, 
x= /(alb—1/4b2) = 2-91 +0-06, 
olf =(1/2b2)4/(4b 4) = 12-82 + 0-42. 


Using the values of a and 6, a curve of | R| against § was computed from (3a)... 
This is shown in the Figure by the experimental points. The Brewster angle, | 
obtained by differentiation of (3 a), to give the minimum value of | R| is found tobe 
79°, as indicated in the curve. The experimental point at 87-5° incidence is in 
good agreement with the theoretical curve, which indicates that there is possibly * 
little difference between the constants governing reflection from water at 8-77 
and 6:35 mm. 

Saxton (1947) gives theoretical curves of the variation of absorption coefficient - 
and refractive index of water with wavelength and temperature. Comparing the * 
above experimental results with these curves it will be seen that the agreement Is. }, 
quite good when the temperature difference is taken into consideration. 
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Coupling of the Ordinary and Extraordinary Rays in 
the Ionosphere 
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BSTRACT. In this paper an approximate solution to the wave equation is given for 
opagation in an ionosphere in which the gradient of the density N is in the vertical, 2, 
irection only, and in which account is taken of the earth’s magnetic field. It corresponds. 
actly to the ray theory and expresses a quantity Z, which is the z derivative of the phase 
ction S, by a quartic equation. Z can be represented as a function of ¢ (which is pro- 
ortional to N) on a four-sheeted Riemann surface, and the branch points are studied for 
e case of vertical incidence for which Z becomes the refractive index. By considering 
e branch points in the complex ¢ plane, the amount of the coupling between the ordinary 
d the isolated extraordinary branches of the (Z, ¢) curves can be expressed as a function 
f the obliquity of the magnetic field. The triple splitting of rays reflected from the 
nosphere, observable where the field is nearly vertical, can thus be explained, and the 
; is substantiated by the observation that the polarizations of the echoes on the (P’, f)- 
cords are ordinary, ordinary and extraordinary in order of increasing critical frequencies, 
: given by the branches of the (Z, ¢) curves. 


HE problem of how a wave travels in the ionosphere is complicated by the 
effect of the earth’s magnetic field, but a very good approximate theory 
can be obtained. It is the more accurate the more uniform the layer is, 
ind it corresponds exactly to the ray theory. 
By the use of Maxwell’s equation we get the original wave equation in the 
vector form 


Vere ed a ; ee LT \e Po vi, A(HE) wy [HE] 

ge Eee al eo 

where vi2=v2(1 +7x), in which vp is the usual critical frequency of the medium, 
4 =v,/2rv the absorption term, », the collisional frequency, v the wave frequency 
y’=v(1+ix), and vy=eH/27m, and is the resonant frequency of the electron in 


the magnetic field H. 
We now assume that the solution of the wave equation can be represented in 


the form of Jeffrey’s approximate solution, 
ant 

~ (dS/dz) 
where E,, E,, E, are the three components of £, and S =coast. is the phase surface, 
and make the approximation that the V2S term can be neglected in comparison 
with the (dS/dx)? terms, which implies that the gradient of the refractive index 


u is everywhere small. 
On reduction, we obtain three equations that can be represented in matrix 


form by 


E Ey 2, 3exp (27uS —2nivt), 


Lm? — 2? +0 hoy X31 
419 1— 27 + G55 mL + Aso x E=0. ac acetae..© (1) 
X13 mZ +- Aog 1 ae m2 + Aga 
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Here we have assumed that the gradient of density in the ionosphere is purely, 
in the vertical z direction, so that we make A(dS/dx)=/ and X(dS/dy) =m 
where / and m give the initial direction cosines of the wave, and A(dS/dz) = a 
011+... are given by 


Oy = Ep + 19 = €pyp2+Inps a1 = €P1P2—IP3 
Moo = Eps + 13 = €1P3 + pr X31 = €p1p3 — MP2 
ogg = fp3 + dog = paps +inpy ago = paps —Inpy 
in which P1,2,3= Ay» 3/H, 
where H,, H,, H; are the components of H, and 
en Vo" Vin. ae Ve Vi (=~ vo" ’ 
py’? — ye, vy? : py’? —y2, yp! (2p? 


From the matrix we obtain an equation for Z which is biquadratic. We can 
thus represent Z on a four-sheeted Riemann surface, and this is sufficiently 
accurate except near the branch, or reflection points. | 

We can express the equation in terms of ¢ and 7, where ¢ is now defined by. 
C=v?/y2, This is the previous ¢ with vy neglected, and +=vy/v, which may be 
resolved into the components 7,,7,, 7;- 

If we know ¢ and the components of 7, we can express Z as a function of it 
position on the Riemann surface (Forsyth 1900). 

We presume ¢ and the components of 7 to be known, because we canno 
attack a problem without knowing the density N in the ionosphere, so that the 
topography of the surface is known and does not depend upon what we represen 
on it. All the branch points and reflection points are known and depend on ¢ 

The behaviour of propagation can be represented as the behaviour of Z, and 
the four-sheeted Riemann surface gives everything that we require. Thus the 
transmission of electric waves in an ionized medium with a magnetic field presen 
can be expressed most suitably in terms of this four-sheeted biquadratic Riemanrm 
surface. This, though not perfectly accurate, is very nearly so, except in th 
neighbourhood of the branch points. 

For the particular case of vertical incidence, the branch points can be determinec 
from the formula 

zap (#8)? = Deo Debate lO) a 
A Raia: ree 


where y=1-—7? and 7,=0. 
Jeffrey’s solution, which gives the waves on the earth, depends upon S, and 
S is given by the above formula, derived from matrix (1). : 
Jeffrey’s solution is : 


A 2nt 
E= asene? (| S,-  2rivt ), oR 


which depends almost entirely on the phase S. 

For vertical transmission, dS/dz can be expressed in the above form. It is 
biquadratic, and represents the variation of the phase of the wave. This formula 
gives the refractive index and all the branch points of S, which are either i 
or infinities of the above expression. 
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We find that the group velocity is zero at all these points, which means that a 
ulse will come up to the branch points, stop, and then go down again. ‘These 
ranch points are therefore the reflection points. We give an example in 
igure 1, which shows the refractive index A(dS/dz) as a function of the quantity 


2, 


, which is the ratio of the frequencies vj/v?, where vj = Ne*/7me?. 


wi Rallies imaginary ~~ = 772 
+ =0-300 Tz = 0-282 Tye = 0°105 


Figure 1. Refractive index as a function of £(« N), showing reflection at the branch points. 


We give here the four branch points in terms of v: 
1p e=ve(1—7) Extraordinary 
2, vet ={(1—72)/(1-22)}!% 
3y ¥, = Vp 
4, y, =voj/(1+7)* 


approximately, since the 7’s are functions of 14, v2 etc. 
More explicitly, they can be expressed in the form 


2 pe} 
a Ash = |p 


1 
1 
5 ko 
a. Vea 0) 
2 2 
+ Vit vgvy = VO: 


* Brought to notice by Rai (1937). 
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The branch or reflection points are shown in Figure 2. This gives the nature 
of the four-sheeted Riemann surface. The branch points correspond to t 
normal ordinary and extraordinary rays. +) 

One thing is clear. The ordinary and extraordinary rays are on different 
sheets of the Riemann surface. When the transmission is vertical and the magneti | 
field is completely vertical, the sheets are in two unconnected pairs, and there isi 
no means of passing from an ordinary sheet to an extraordinary one or vice versa. 
There is no coupling between the rays. This is discussed in fuller detail later 
but this outline gives an idea of the nature of the coupling between the two. 

We can now go on to consider the Riemann surface method of representation 
The results contained in the following pages give us an understanding of the 
coupling process. No criterion of the relative probability of appearance or the 
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Figure 2, Refractive index as a function of ¢(« N), showing branch points in the complex ¢ plane¢ 


degree of intensity of the rays or pulses of these critical frequencies has ever been 
given, and the following method of analysis gives definite conclusions with regardk 
to this point. ‘The original analysis, giving critical frequencies alone, is like the: 
older phase integral quantum theory, in which the various spectral lines were givenr 
(corresponding to the critical frequencies) without any indication of their intensitiess 

We have represented the waves on a four-sheeted Riemann surface on whichi 
the branch (or reflection) points will be seen. It has been pointed out that the 
phase and group time increase continually up to and through the branch poiats: 
and follow on the downward branch to ¢=0, giving the complete phase and group; 
delay of a wave or impulse travelling along the upward branch (dS/dz) positive’ 
to the reflection point (d.S/dz)=0 (generally but not necessarily) and down the 
lower branch (dS/dz) negative to the starting point. 

The conditions of reflection can therefore be best studied by examining the: 
branch points and branch lines of the four-sheeted Riemann surface on whic 
(dS/dz) is most suitably represented. 
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This surface, for the case discussed above, is shown in Figure 2, where the 
rface represents the complex { plane. ‘The study of expression (2) shows that 
anch points occur at 


F=1—- (one), ¢=1 (one), {=1+(ér?/2r,) (two), 
C=1 —(t7r2/27,) (two), €=1+7 (one), £= co(unless (= 27), 
and £=y/(1—72) (one), 


hen the refractive index becomes infinite. 

If £ negative is cut off by a perfectly conducting plane, there are two further 
ranch points at =0 representing the complete reflection of the ordinary and 
xtraordinary rays at the surface. The branch lines ending at the branch points 

ould be noted. 

Any path which crosses a branch line goes from one sheet, say 1, to another 
heet, say 2. The branch line ends at any branch point where the two are 

terchanged. 

The analysis of the four-valued expression for 4 ic this manner reveals signi- 
cant features. When the incidence is vertical and the earth’s magnetic field is 
so vertical, then there are two independent solutions on two independent 
airs of sheets, i.e. the 1-2 pair and the 3-4 pair, and there is no interchange 
etween the pairs. 

: When, however, the magnetic field is at all oblique, both the longitudinal and 
the transverse components of H are present. All four sheets are then linked 
together. There are branch points where 1 and 3 and where 2 and 4 are inter- 
changed, and it is always possible to choose a path which goes from any one sheet 
to any other. The ordinary and extraordinary rays are then coupled together, 
as well as the ordinary and isolated rays, although usually to a very slight extent. 

This analysis may be used for estimating or computing the intensities of the 
echoes corresponding to the four critical frequencies already given, and for 
calculating the coupling between the ordinary and extraordinary rays. 

In the first place it should be noted that it will be very improbable that the ray 
with critical frequency v2 (pointed out by Rai 1937) will ever be obtained with 
observable intensity, for the path of the extraordinary echo which is reflected at 
£=y/(1—72) (see Figure 2), corresponding to the critical frequency v’, has to 
traverse the region from =1—7 to C=y/(1 —7%), where the refractive index is 
large and entirely imaginary, before being reflected. 

There is, however, a much more likely case, in which, instead of the normal 
ordinary and extraordinary doublet, which separates as a split echo near the 
critical frequency, there is a triplet ordinary, ordinary, extraordinary: 0., 0., €X. 
(see Figures 2, 3 and 4). 

The three critical frequencies here shown are given by 

¢=(1+7), C=1 and {=(1—7). 
0. O. ex 
Since € = v2/v, where v is the critical frequency, the three critical frequencies 
are 
vo/(1 +7)3, Vo and vy/(1 —7)? 
0. O. ex: 
where v2= Ne?/7mc?. 
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The last two are the normal critical frequencies and the first an extra one. 
This split into three branches was first observed by Harang (1936), who 
showed the three branches (a), (b) and (c) each separated by a frequency interval | 
v_/2 appropriate to the critical frequencies given above. oP | 
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Figure 3. Similar to Figure 1, showing the corresponding (h’, v) curves with the 
production of a triplet. 


h 


CBA 
aoe oe 6 
A 8B C 
y Echo Pattern corresponding to 
(a) 


Frequency shown in (a) 
a (b) 


Figure 4. (h’, v) curve, showing the triplet with the corresponding echo pattern and polarizations. |, 


The mechanism of this is illustrated in Figures 1, 2 and 3. Turning first to | 
Figure 1, it will be seen that if we could bridge the gap from the ordinary ray (1) | 
to the isolated ray (3) at C, and then move along the 3-4 path through the apex 
at D to C’, bridging the gap again at C’ on to the downward ordinary branch 2, 
we should form a path the critical frequency of which would be Vo/(1 +7), 
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rresponding to the added ordinary ray observed by Harang. ‘This ray was also 
seryed by the author oa the cathode-ray oscillograph in 1933, the order of delay 
ing ex., 0., 0. (Eckersley 1933, also Meek 1948, Newstead 1948). 

The polarization of this extra ray, determined as it is by the conditions when 
aving the ionosphere on the downward ordinary path, will obviously correspond 
that of the ordinary ray. 

That the physics of such a transmission is possible will be seen by proceeding 
the limit where 7,->0 and H, the magnetic field, is purely vertical. When 
, is very small, but H not quite vertical, the analysis shows that the isolated 
ranch 3-4 very aearly touches the ordinary ray branch at CC’, and in the limit 
hen 7,0, the two actually join, and the paths CAC’ and Coo disappear. 
e are then left with the two simple rays with critical values of ¢, viz. 1—7 and 
47. These are the extreme rays of the triplet, and have a frequency difference 
y instead of the usually observed difference vq/2. 

It is hardly conceivable, physically, that this process should occur discon- 

inuously at r,=0, and we should expect, on physical grounds, to find all three 
yranches in the neighbourhood of 7,=0 for the reason that the isolated branch 
34 is closely coupled to the ordinary branch 1-2 when 7, is small. 
We can express these ideas in a more exact form by considering the Riemann 
surface again. We have to find the total attenuation on integrating this along a 
path which, starting in sheet 1, passes continuously to sheet 3 round the 1-3 
branch point, again passes continuously from sheet 3 to sheet 4 round the 3-4+ 
branch point at €=1 +7, and then from sheet 4 to sheet 2 round the 4-2 branch 
point and back to earth. In sheet 2, as in sheet 1, N, the number of electrons, 
is zero at C=0. 

This represents, na moreexact mathemetical form, the physical path previously 
discussed. ‘Thus it passes, in sheet 1, as far as C=1; then it goes from 1 to 
1+7(72/27,), where, ¢ being complex, there is an attenuation. After passing 
round the branch poiat 1-3, it reaches sheet 3, and ¢ decreases on this sheet from 
1 +2(72/27,) to the neighbourhood of ¢=1, again making a contribution to the 
attenuation. We then integrate along the path from t=1 to €=1+7 in sheet 3, 
where the path is along the real axis and the integral is real, so that there is no 
contribution to the attenuation. At this point, 1+7, it goes from sheet 3 to 
sheet 4. 

Finally, we require to pass from sheet 4 to sheet 2 to continue the path 
down to earth, and to do so we have again to introduce attenuation by rounding 
the branch poiat at C=1 —i(r2/27,) where ¢ and the integrand are again both 
complex. We then reach branch 2, the ordinary ray, at €=1 and proceed, without 
attenuation, along this path to earth again at C= 0) 

The imaginary part of pu is definitely specified at each point of the path, so 
that the attenuation can be uniquely determined. 

It will be observed that the path and integrand are only complex in the regions 
between (=1 and €=1 +1(72/27,) and between (=1 and €=1—1(7;/2r,), and that 
this path tends to zero as +, 0. Therefore the total attenuation of this spurious 
ordinary echo decreases as T,—> 0. ‘This was also obvious in the previous more 
definitely physical treatments. 

Another point must be considered: The integral should be taken with respect 
to 2, i.e. the distance, but as expressed on the Riemann surface, p, or A(dS/dz), 


is a function of ¢ only. 
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The attenuation or imaginary part, 
dS ei 
[rgeaz= [A dz, | 


can be expressed in the form [f(¢) d(dz/d0). = | 

Now the first part of this integral is a quantity only depending on 7,, 7,, 7, 
not on the distribution of density in the layer. The whole integral is small when) 
the gradient of density in the layer is large, and vice versa. 


vy 


Therefore the attenuation equals [ KO fea , and is likely to be small when 


d{/dz, the ionic gradient, is large, and vice versa. 

It follows that the conditions in which the triplet form with critical frequencies 
vo/(1—7)?, v9, v9/(1 +7)? is likely to appear are: (i) in northern latitudes where 
T, 1s small, and (ii) in regions where the gradient of N with respect to ¢ is large. 
It is significant that this triplet form has only been observed in high latitudes. 

Cases which may well represent this type of transmission were noted in 1933 
(Eckersley 1933). On such occasions well defined triplets were observed on 
frequencies in the neighbourhood of the critical frequency. ‘The polarization 
observed in the cathode-ray indicator was invariably as shown in Figure 4. 

A and B are the normal extraordinary and ordinary pair, the extra echo C! 


being of ordinary polarization. 4 
These three correspond to the branches with critical frequencies : 
A B S 
vo/(1 —7)3, Vo and vo/(1+7)* (the extra echo). 
ex, 0. O. 


We may use the curves in Figure 1 in the way just described to calculate the © 
coupling coefficient between the ordinary and extraordinary rays. Thus, 
considering an up-going ray in the ordinary 1-2 sheet, we may bridge the gap, as 
before, to the isolated ray in the 3-4 sheet. We then pass to the branch point at 


C=y/1—-7,)!, dS/dz= x 


and back along the imaginary branch in the 3-4 sheet to D, and finish along the © 
real part of the extraordinary branch which goes down to earth again. It is clear | 
that the total attenuation is large and, therefore, that the coupling will be small. — 
The contribution to the attenuation is especially large between Sand D. Thesame | 
remarks as regards the gradient of C with respect to s apply equally well here. Thus | 
if d¢/dz is large enough, the total attenuation along the path may be made as small 
as we please. 

The truth of this statement may readily be seen by proceeding to the limit 
where d¢/dz becomes so large at one point that we may consider this to be the i 
surface of separation. This case may easily be worked out according to the full © 
wave theory. Without going into detail, it may be said that in medium 1 (see __ 
Figure 5) there may be two characteristic waves, 0., and ex.,. Let the polarization | 
ratio X/Y be given by P,,,, (pect | 

Therefore, in medium 2, we shall have Veet agitate 

Suppose a wave of the type o. is sent through the first medium. Since 
P,2 and P.... are not the same as P,, and P..,, and since there is an appreciable 


| 
| 
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component of the electric force, the reflected wave will not be only of the o. type, 
ut will consist of o. and ex. ia suitable proportions to satisfy the boundary 
nditions. 
If the path in medium 1 is sufficiently long, and if the incoming signal consists 
a short wave impulse, the reflected wave will consist of an 0. and an ex. impulse 
hich, after travelling sufficiently far, will be separated into two pulses. 
Thus, by sending into the medium an o., impulse, we get a pair, Ry, and 


ex reflected or refracted back into medium | again. 


: Transmitted 
Waves 


0 
Reflected 
Waves 
e 


5 Direction 


Figure 5. Illustrating the full wave solution at a surface of separation between ¢; and 2p. 


Thus there is obviously a coupling between the o., and ex.; rays where the 
gradient is so steep that in the limit we may consider it to form a definite surface 
of separation. 

This is an extreme illustration of the effect of coupling already referred to. 

This coupling has, in general, been calculated by the relatively inaccurate 
four-sheeted Riemann surface method because it is very difficult to get a suitable 
solution of a fourth order differential equation. As reckoned by this method, 
the amount of coupling may be a little inaccurate, but where the media on each 
side of the surface of separation are uniform, it is possible to use a completely 
definite wave solution which agrees with the above and shows, as in the last analysis, 
that there is a definite coupling between the o. and ex. waves. Some German 
papers also show this, and I think we may rest assured that this coupling exists. 


RECAPITULATIONS 


1. Triple splitting of rays reflected from the ionosphere has been observed 
in the polar regions where the earth’s magnetic field is nearly vertical (see Figure 4). 
Well defined triplets are almost invariably of the type 0., 0., ex. in terms of 
critical frequency. In terms of delay time they are ex., 0.,0. "They have been 
observed by Harang in Norway, J. M. Meek in Canada and by the author as far 
south as England. They have also been observed by Newstead in Hobart, 
Tasmania, but as far as is known they have never been observed in any equatorial 
region. 

2. These triplets can be explained in terms of the four-sheeted Riemann 
surface, which expresses the propagation in the ionosphere. They are due, as 
shown in Figure 3(6), to the coupling between the isolated and ordinary rays 
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This coupling does not occur in equatorial regions. ‘The Riemann surface theory, 
which is approximate, has been confirmed in certain cases by the wave theo fl 
which is absolutely correct. | 

3. It is in virtue of a coupling between the o. and isolated rays that a triplet is; 
formed. When 7, is zero, and we are consequently at the pole, there is no definite ¢ 
isolated ray, and the two extremes of the triplets are formed with a frequency| 
separation of yy, approximately, aad not the usual v,/2 approximately. 

4. The fact that the observed features of the triplets can be explained by thet 
Riemann surface theory, and can be derived fron it, is a very good proof that 
this theory is correct. These features are: (i) the triplet is of the ex. 0., 0. type; 
(11) it is confined to high positive and negative latitudes. 

Triplets were noted by the writer in England as early as 1933; experimental 
results were given, but the explanation was not published at the time. 
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EEETERS “TOSTHESEDITORN 
The Debye Effect in Electrolytes 


P. Debye (1932) has predicted the existence of an alternating electric potential ® at a4 
point in an electrolyte through which a sound wave is passing. His equation is 


pala SPS EGR 
«°° 3p; 6;7/p;) \[1+G4nl/Da)*} 


where w is the angular frequency, g the velocity of sound in solution, ao the particle velocity , 
amplitude, / the conductivity of the solution, D the dielectric constant of the solvent, My 3 
the mass of a hydrogen ion, « the electronic charge, p; the relative number per unit volume 2 
of ions of valency ¢;, gram ionic mass M; and frictional coefficient p;. "The deduction of f 
this result has been followed by mathematical work by which corrections for various smaller t 
effects have been proposed. 

For frequencies less than tens of megacycles per second the last factor in the equation 1 
is of the order of unity and frequency becomes unimportant. For a binary salt with 1 
univalent ions Debye’s result then reduces to 


M,/pi—M,/p2 
1/pi+1/pe 


With dilute solution in a given solvent, electrolysis experiments give values for ionic mobilities 
differing little from ion to ion or with concentration. Hence for the purpose of determining : 
the order of magnitude of the effect, we may put Pips so obtaining a result roughly / 
proportional to (M,—M,). ‘The order of ® should therefore be about 1 microvolt when } 
a 1s 1cm/sec.: a sound intensity of about 1/100 watt/cm®. This led Debye to express 5 
confidence in the early detection of the effect. In fact, no experimental evidence for it in) 
electrolytes was found until the recent work of Yeager, Bugosh, Hovorka and McCarthy / 
(1949), although an effect of about one hundredth of the estimated order of magnitude was.} 
recorded by Rutgers (1946) for colloidal solutions. | 


@=1:4x10-‘ay volts. 
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For purposes of calculation we should use for M, and M/, the effective masses of the ions, 
cluding the mass of any solvation sheath. ‘The use of the available conflicting results for 
lvation numbers for calculation of the relative values of ® for the five salts used by the 
ove workers shows no correlation with their experimental results. In fact these results 
ry so little between one salt and another in water as to suggest that the solvent is largely 
sponsible for the effective masses of the ions. That is, there is either a large solvation. 
eath masking the effect of the ions proper, or the equation does not take into account all 
e factors involved. Available results from which conclusions may be drawn are insufficient 
d work with selected salts and a variety of solvents is obviously very necessary. Such 
sults should serve to obtain quantitative data, not only about the Debye effect, but also 
out the mechanism of solvation. 

A principal experimental difficulty is the possibility of direct electromagnetic coupling 
etween the radiofrequency source and the detecting apparatus. Such coupling would 
mit the sensitivity of the amplifier which could usefully be employed. Yeager, Bugosh, 
ovorka and McCarthy used elaborate screening precautions to detect the effect and their 
ethod thereby made measurement of the sound intensity impracticable. However their 
sults appear to confirm that the effect is smaller than expected from the Debye equation 
nd for reasonable sound intensities is therefore not much above noise level. 

A solution would seem to be the use of a method by which a pulse of sound is sent 
Bough the experimental liquid, the duration of the pulse being such that, with the length 
f sound path used, the source is inoperative during detection of the effect. 

An apparatus of this kind, by which any electromagnetic pick-up is separated in time 

ss the desired signal, is nearing completion in this laboratory. The frequency of 

165 kc/s. has been chosen so that the sound wavelength in the experimental liquid may be 

sufficiently long for the use of probes spaced half a wavelength apart for detection of the 

Dotential variations; at the same time standard intermediate frequency transformers may be 

ised for the amplifier. At this frequency the absorption of sound in the electrolyte will be 

small and a long sound path may consequently be used. This makes it possible to use a 

ong pulse, of the order of one millisecond in duration. The amplifier required for such 

1 pulse need have a bandwidth of only 1,000 c/s. and the noise level will therefore be low. 

The method is essentially a progressive wave method, reflections from the end of the tank 
being avoided by the use of an absorbing device. Since a large tank may be used the sound 
Geld at the detecting probes is calculable and no difficulty is expected in determining the 

intensity there. The pulses are to be displayed on a cathode-ray tube and compared with 
signals from a standard generator. It is hoped to investigate thoroughly the dependence of 
the effect on the various parameters involved. 


Physics Department, A. N. HunrtTER. 


University College, 
Leicester. 
7th October 1949. 
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Electron Optics of the Three-Stage Electron Microscope 


Electron microscopes employing three stages of electronic magnification are now 
becoming more generally used, but as yet no adequate account of their electron optics has. 
been published. Le Poole (1947) gave some ray diagrams, which, while they are accurate,, 
are not sufficiently comprehensive to give a clear picture of the situation. 

In connection with the building of an electron microscope in collaboration with Mr. J. F. 
Brown, and under the direction of Professor G. I. Finch, the conditions of conjugate foci 
between the electron lenses were examined by Gaussian optical methods, which, while they 
do not give a strictly accurate quantitative analysis of the situation, do indicate the approxi- 
mate positions and method of formation of the images. The analysis is reproduced in the 
Table, and Figures 1, 2, 3 and 4 are the ray diagrams which were deduced from the data 
summarized in this Table. 
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Condition of Conjugate Foci between Objective Lens and Intermediate Lens of a 
Three-Stage Electron Microscope. 
u (objective lens)=0-8 cm.; v (intermediate lens) =17 cm.; objective-intermediate 
lens separation=15 cm.; magnification of projector lens = M,=120. 
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Figure 3. Low magnification (a). Figure 4. Low magnification (). 


Figure 1 indicates the electron paths when the microscope is producing a high magnifi- 
cation final image. If from this condition the power of the intermediate lens is decreased 
sufficiently, it is possible to obtain a zero magnification image, and if no objective aperture is: 
present to remove the higher order diffracted rays, a diffraction pattern of the specimen is 
seen on the final screen (Figure 2). If the power of the intermediate lens is further decreased 
and the objective lens refocused, then a low magnification image is produced on the final 
screen, and the optics of this image are best explained in the following way. 

If the intermediate lens is switched off, and the objective lens refocused, a two-stage 
image is formed. If now the intermediate lens is switched on at very low power, then the 
magnification of the image is slightly decreased ; this condition is shown in Figure 3. Upon 
further increasing the power of the intermediate lens, a stage is reached where the objective 
lens no longer produces a real image ; this is the case illustrated by Figure 4. Further 
increases in intermediate lens power produces again the diffraction pattern. 

These theoretical deductions were borne out by practical measurements of focal lengths 
of the lenses of the electron microscope for the various conditions of image formation. It 
was also found that, when a low magnification was required, the image formed by using the 
intermediate lens in the condition corresponding to Figures 3 or 4 showed less distortion 
than that produced when the conditions of Figure 1 obtained. 

I am indebted to Mr. J. F. Brown, with whom the work was carried out, and to Professor 
G. I. Finch for his helpful supervision. 


Applied Physical Chemistry Laboratory, Gaba CHrArnicr= 
Imperial College, London S.W. Hee 
20th September 1949. 


Le Pootg, J. B., 1947, Philips Tech. Rev., 9, 33: 
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Metal Rectifiers, by H. K. Hentscu. Pp. xi+156. First Edition. (Oxford: 
University Press, 1949.) 15s. 


This book belongs to the series of monographs on Physics and Chemistry of Materials. 
I always find it difficult to define exactly the cross section of physicists or engineers to whom) 
‘such a monograph is addressed. In general the practical expert will find little new in the, 
-chapters on manufacture, while he might find the chapters dealing with theory hard going 
because they are, of necessity, rather condensed. On the other hand the theoretical 
physicist often finds that the description of the practical methods is somewhat incomplete, 
It appears to me that in this volume the author has managed to find a happy combination 
of describing the essentials of practical manufacturing methods without getting lost in details, 
and of presenting the theory in such a way that it is not too difficult to follow for the 
practical engineer, though I myself would have liked to have seen the theoretical Chapter V 
somewhat amplified and some of the figures in it described in more detail. On the other 
hand, a full bibliography enables the interested reader to study the subject in greater detail.| 

After a short history of rectifier development the manufacture of the most important types 
-of rectifiers is described. These comprise cuprous oxide, selenium, sulphide and germanium 
‘rectifiers. 

In the third chapter characteristic rectifier properties are described. Starting with the 
rather complex voltage—current characteristics of the various types of rectifier, D.c. resistance 
and incremental or A.c. resistance are described as well as the phenomena of hysteresis, creep, } 
ageing and the effects of temperature. The self capacitance represents an A.c. shunt and az 
number of equivalent networks are given together with an account of the variations of the self 
capacitance with the conditions of measurement such as direct voltage, time, etc. In view 
of the importance of “‘ noise’ in modern circuitry, the question of noise would appear to 
deserve more than one short paragraph at the end of the chapter. 

Chapter IV deals with measurements of rectifier properties and the pitfalls that may be 
encountered in these. 

The longest chapter in the book is, quite rightly, the one dealing with modern theories of { 
rectification. After a discussion of the structure of semiconductors and the mechanism of } 
electrical conduction in solids, modern concepts of insulators, intrinsic semiconductors and 
impurity semiconductors are given. 

Next, the theories of the chemical and physical barrier layers are explained and the : 
theoretical and experimental results compared. Chapter VI deals shortly with earlier 
theories of rectification. 

Two short chapters on rectifier operation and future developments conclude this highly ; i 
recommendable book. As far as the bibliography is concerned, replacement of the simple 
asterisks by numbers or by an indication of the corresponding page in the text would | 
facilitate finding one’s way back into the text when scanning through the list of references. . 
Seven of the references mentioned in the text have not received an asterisk, and the one » 
mentioned on page 129 has the wrong number. H. G. LUBSZYNSKI. - 
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Heat Transport in Superconductors, by K. MENDELSSOHN and J. L. OLSEN. 


ABSTRACT. The heat conduction processes in superconductors have been discussed 
on the basis of the analogy with liquid helium II. It is suggested that under certain con- 
ditions the superconductive metal can exhibit a type of heat transport corresponding to the 
heat flow associated with the fountain effect. The heat conductivity of a number of pure 
metals and alloys has been measured in the normal and in the superconductive state and the 
results have been analysed with reference to the hypothesis mentioned above. In addition 
to the change of heat conduction with temperature, the magnetic hysteresis of the heat 
conduction has been investigated. The use of these phenomena as make-and-break 
thermal contacts at very low temperatures has been suggested. 


The Magnetization of Superconducting Plates in Transverse Magnetic Fields, by 
E. R. ANDREW and J. M. Lock. 


ABSTRACT. ‘The magnetization curves of thin superconducting tin plates were measured 
in transverse magnetic fields in order to investigate the nature of the intermediate state in 
such specimens. The curves showed peak magnetizations much higher than H,,/47, indi- 
cating the presence of fields much greater than the critical field at the edge of the plates, 
This effect is qualitatively similar to that predicted by Landau’s theory of the intermediate 
state, but quantitative agreement is not obtained. Resistance measurements on thin strips 
of tin in transverse fields showed that resistance only reappears for fields considerably 
greater than those required to start the destruction of superconductivity. 
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A Modified Perturbation Procedure for a Problem in Paramagnetism, by M. H. LJ 
PRYCE. Hi 


ABSTRACT. A modified perturbation technique is described for problems in whie 
second-order effects are comparable in magnitude with first-order effects, where orthode 
methods break down. It is applied to the energy levels of paramagnetic ions in a crysté 
giving an effective Hamiltonian in which the Stark splitting, the anomalous g-value and tl 
temperature-independent paramagnetism are clearly exhibited. 


The Paramagnetic Resonance from Nickel Fluosilicate, by R. P. PENROSE a 
K. W. H. STEVENS. 


ABSTRACT. The paramagnetic resonance from the nickel ion in nickel fluosilicate “at 
been investigated at a number of frequencies and temperatures. It has previously bee nt 
suggested that the ion is in an electric field of trigonal symmetry. The results obtained are § 
in good agreement with this hypothesis. It is also found that the magnitude of the splits 
of the levels varies with temperature. : 


The Theory of Magnetic Rcsonance-Line Widths in Crystals, by M. H. L. PRYCE 
and K. W. H. STEVENS. e 


ABSTRACT. A theory is developed for calculating the mean displacements and mean) 
square widths of resonant absorption lines in crystals for which the spin-lattice relaxation: 
time is long compared with the spin-spin relaxation time. ‘The theory is applied to a¢g 
number of different cases ; in particular, the effects of nuclear hyperfine structures and 
exchange forces in ionic resonances are discussed at length. It is shown that the theory ; 
holds, provided that the temperature is high compared with the Curie temperature. 


On the Capture of Particles into Synchrotron Orbits, by 'T. R. KatsEr. 


ABSTRACT. A theory is developed for the capture of particles into synchrotron orbits: 
when the radio-frequency accelerating voltage rises from zero to its maximum value in a 
finite time. It is shown that if this time of rise is sufficiently long, all particles originally 
occupying a band of width 1/1/2 times the maximum width of the final stable region for 
synchrotron phase oscillations will be captured, irrespective of the initial phase at which they ; 
enter the accelerating gap. ‘This conclusion would seem to be in agreement with thet 
observed characteristics of existing electron synchrotrons. 


Experiments on Electron Capture and Phase Stability in a 14 mer. Synchrotron, | 
by TI. R. Kaiser and J. L. Tuck: 


ABSTRACT. Experiments were made by interrupting the radio-frequency accelerating 3 
voltage for intervals during the acceleration in a synchrotron, and observing the effects on ¥ 
the output. The process of capture is found to be more efficient than would. be expected f 
for an instantaneously developed radio-frequency accelerating voltage. The synchrotron 
acceleration may be interrupted for a short period without losing more particles than would | 
be accounted for by contraction of the beam radius during the interruption, so that for 
sufficiently short periods of interruption the loss is negligible. "This demonstrates that the 7 
phase at which electrons enter the resonator during capture is of no consequence under the ? 
conditions prevailing. ‘The results are in quantitative agreement with a theory put forward ! 
by Kaiser which takes into account the finite rate of rise of resonator voltage, and with the 
general theory of synchrotron stability. 
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